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ABSTRACT 


The Frequency Planning Problem, for large communication 
netv>c>rhs (like Indian Railways Microwave network) is considered 
where the interference-free operation of the network is managed 
by a one-time assignment of discrete frequencies to the various 
links. The problem can be formulated as a nonlinear integer 
program. However, the nonlinear integer program is difficult to 
solve for a reasonable size problem. Alternatively, the problem 
can be modelled as a Graph Coloring Problem with some additional 
constraints due to adjacent— channel and intermodulation 
interference mechanisms. This approach was found feasible and is 
followed. The graph being reasonably large and the Graph 
Coloring Problem being NP-Complete, some heuristic algorithms are 
sought to solve this problem. 

Sequential algorithms for Graph Coloring are considered, 
with modifications to take care of the additional constraints due 
to adjacent-channel interference. A scheme of partitioning the 
graph into smaller sub-graphs, which can individually be colored 
quite easily and efficiently, is suggested, since interference 
mechanisms have a local effect. Next, the individually colored 
sub-graphs are put back together following a systematic method of 
reconciling the conflicting colors to . get back the original 
graph. The frequency assignment corresponding to the graph 
coloring is evaluated for i nt er-modu lat i on interference. 
Critical cases of interference are identified link-wise and it is 
left to the user to modify the frequency assignments iteratively 


and eliminate these cases. 
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CHAPTER I 


INTRODUCTION 


In this thesis, an attempt shall be made to provide a 
methodology for fr&Quency planning for large cormcunication 
networks, keeping in view the requirements of a typical network 
such as the Indian Railways Microwave communication network. In 
order to understand the nature of the problem, a brief 
introduction is given in this chapter, to the subjects of 
communication networks, frequency planning and 
interference-caused constraints on frequency planning, followed 
by a survey of the literature dealing with frequency planning in 
particular. Next, the frequency planning problem is explained in 
brief. An overview of the thesis and, its organization is 
provided in the final section. 

1*1 Communicat-ion Nat works and Frequency Planning 

In today's era of high technology, electronic 
communication technology is seen at the forefront. Communication 
networks span entire nations and most of the world, relaying a 
broad spectrum of signals such as data, voice and picture from 
one place to another. The single most precious resource used up 
by these networks is, perhaps, the electro-magnetic spectrum; 
its useful range being determined by the technology put to use. 



For example, microwave networks will have a useful spectrum range 
from 5 GHz to 50 GHz, while domestic radio broadcasts use a much 
lower range of frequencies. Frequency planning is directed 
towards assigning different segments of the useful range of 
e 1 ectro-magnet i c spectrum to different links of the commun i cat i on 
network with the objective of ensuring a reliable, 
interference-free performance of the network. Frequency planning 
is therefore a vital exercise towards the efficient utilization 
of electro-magnetic spectrum by a communication network. 

A communication network can, very s impl i st i cal ly , be 
described as an undirected graph G<N,A) where N is the set of 
nodes which are nothing but the sites for transceivers 

< transmi tter-cum-rece i ver) and A is the set of arcs joining two 
nodes in N, if there is a direct communication link between the 
two nodes. (See Figure 1.1) The actual operation of a link in a 
communication network is described in Appendix I. 

In actual practice, the sites may be the metropolitan 
cities of India and the links may be the satellite links relaying 
telephone and television signals between the cities. A more 
elaborate instance of a communication network, frequency planning 
of which provides a background to this thesis, is that of the 
Indian Railways (IR). which uses the microwave (MW) range of the 
e 1 e c t ro-magne t i c spectrum. The frequency planning problem for 
such a network is explained later in this chapter (Section 4). 






FIGURE 1.1 Representation of a communication network by a graph G 
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Frequency planning for such a network implies, firstly, 
the creation and maintenance of a database of the transceiver 
sites and the communication links including all their relevant 
characteristics such as geographical location, transmitter power, 
and receiver gain. Secondly, frequency channels are assigned to 
links such that compatibility in terms of electro-magnetic 
interference (EMI> is likely to be satisfied. Next, calculations 
are made for the compatibility of each link, in terms of EMI, 
with the already planned network. If compatibility is not 
satisfied for a link then the frequency assignment to this link 
IS modified. When compatibility gets satisfied, the frequency 
assignment to this link is documented and the planned network is 
updated. An extended overview of the frequency planning activity 
in Indian Railways, for their MW network, is given in Appendix 
I I . 

In addition to the above, frequency plans are subject 
to obeying the international regulations on the use of the 
electro-magnetic spectrum such as those implemented under the 
International Telecommunication Union (ITU) or CCIR. Moreover, 
at the national level, co-ordination is required between the 
various users using the same range of the spectrum. 



1.2 Constraints on Frequency Planning due to Interference 
Mechanisms 

Re-stating the objective of a frequency plan, it is to 
ensure that EMI is kept within tolerable limits for the operation 
of each link of the communication network. In this context. it 
will be useful to identify the main sources of EMI. 

According to the classifications described, for 
example, by Zoellner and Beall Cl], constraints imposed by EMI in 
the frequency plan of any communication network may be classified 
into three groups-, namely, 

<i> Co-chanriBl Constraints Certain pairs of transmitting 
links cannot be assigned the same frequencies due to 
the locations of antennae, power of signals, and other 
characteristics. 

(ii) Adjacent -channel Constraints Certain pairs of 

transmitting links should be assigned frequencies that 
are separated by a minimum amount in order to avoid 
i nterf erence . 

(iii) Co-silte Constraints Particular subsets of 

transmission frequencies should not be assigned to 
convnun 1 cat 1 on links which are located in close 
proximity in order to avoid potential interference such 
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as i nternnodulat i on (explained in Chapter II). 

Appendix I provides a short introduction to the above 
interference mechanisms which give rise to the above constraints, 
how they arise and how their magnitude can be estimated by 
various formulae. 

Technological means exist to reduce the undesirable EMI 
in both transmission as well as reception. Band pass filters do 
not allow signals outside of a specified frequency band to pass 
through. Antennae are of di rected-type so that s i gnal -spread is 
less in other directions. Polarization of signal waves is also a 
useful option to reduce the interference effects. Some of the 
above techniques receive a greater attention in Appendix I. 

Given a hardware setup for operating the communication 
network, including all the technological means quoted above, the 
interference levels depend upon two factors 

(i) Geographic a.1 Locatior. o/ t?\e Links If an antenna of 

one link "sees” the antenna of another link and the 
length of the signal path between the two antennae is 
significantly low then the two links can potentially 
interfere with each other’s operation. 

<ii> The Frecfxiency Plan for the Netxoorh If two links are 
located such that they can potentially interfere with 
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each other then larger the separation between the 
frequencies of the two links, lesser is the actual 
amount of interference. So a good frequency plan takes 
care in assigning frequency channels to such links. 

There may be other spurious sources of interference in 
the environment and the measure of the actual EMI may be 
statistics like the average number of interference cases per 
month, mean time of interference per month or any other suitable 
measure for the application at hand. In general, the 
interference function will be highly nonlinear both in the time 
domain as well as the frequency domain. 


1-3 Litera-ture Survey in the area of Frequency Planning 

Mathematically, co-channel and ad j acent -channe 1 
constraints are relatively easy to describe in comparision to 
co-site constraints. Frequency planning scenarios with these 
constraints have drawn some attention. Several approaches have 
been proposed for these type of problem scenarios. 

Freeman C2], in a Dynamic Prograrmiirc^ formulation, 
suggested building the optimum frequency plan in a stage wise 
manner by minimizing a certain Interference Potential Index. 
Some of the indices he suggests are Oeometric Mean Index 
(geometric mean of Interf erence-to-No i se ratios over all links). 
Mean Interference Poxoer Index (simple arithmetic average). 
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Maximum Interference Power Index (maximum over all the links) and 
Interference Cases Index (number of interference cases in the 
environment). The notion of Expected Valxie Index of interference 
is also suggested for a more rational approach. 

A Mutual Interference Chxxrt method due to Sachs C5] 
gives an analytical concept for doing compatibility calculations 
before making frequency assignments. 

Zoe liner C4] has proposed a Simulation Game Theoretic 
approach to evolve strategies for frequency assignment. A 
Multiuariate Search approach to frequency assignment using 
Nonlinear Programming has been given by Massaro C5]. Cameron C6] 
has treated the problem as a Graph Coloring Problem, as well as a 
Set Covering Problem. 

Thuve C7] has also recast the frequency planning 
problem as a Set Partitioning Problem with the objective of 
minimizing the Total System Interference . He has however used a 
Link Interference Function which is a combined measure of the 
transmitter and receiver interference at the link terminals. 
This function is, in general, spiky and highly non 1 i near , but the 
solution methodology is independent of it's nature. The 
co-channel and the adjacent-channel constraints are handled as 
frequency separation constraints. 
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More recently, the co-site constraints arising due to 
i ntermodu 1 at 1 on interference have received some attention in the 
frequency planning problem. Mathur et al C8] derive an Integer 
Nonlinear Programming model to analyze the extent of 
intermodulation interference at a site in a communication 
network. They also extend the algorithms developed using this 
model to generate Intermodulation-Free subsets of a given set of 
frequencies which is useful in making compatible frequency 
assignments. Based on the same model, Mathur et al [9] have also 
described an Interactive Computer Software System for the 
frequency planning problem. 

Box CIO] describes a heuristic frequency assignment 
technique in which the links that prove themselves to be 
"difficult” in respect of interference-free assignment, rise 
rapidly towards the top of the list of links. The iterative 
method can be easily automated- However, for cases where no 
frequency assignment is possible with the existing resources, 
this method may not terminate quickly. 

A survey of the applications of quantitative techniques 
in Radio Communication Interference problems is provided by 
Morito et al Cll]. 

In a related area, Hafez et al [IZ] examine the effect 
of user's geographical distribution on the utilization of radio 
spectrum. In general, the remaining spectrum available for 
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further assignments at any time is usually less than the unused 
portion of the spectrum. Using appropriate measures of Spectrum 
Utilization Efficiency, frequency plans and strategies for 
frequency assignments can be evaluated for their efficacy. In 
this context, the authors present a new measure of true spectrum 
utilization and develop a model to calculate it. 

1.4 The Frequency Planning Problem for Indian Railways Microwave 
Communication Network 

The IR-MW network has about 100 sites and 500 links 
spanning the entire Indian nation. A brief description of 
frequency planning activity for Indian Railways MW network is 
given in Appendix II. Each link in the network operates on a 
discrete frequency channel. As the number of frequency channels 
available for use is limited, they have to be shared between the 
large number of links. A signal transmitted along a link may 
stray into the receivers of a neighbouring link and cause 
interference problems there. The magnitude of interference is 
significant whenever both the links are operating on the same 
frequency channel or adjacent frequency channels. The frequency 
planning problem (FPP) is to allocate the frequency channels to 
the links such that there are no significant interference 
problems for any of the links. 
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1.5 Overview of t-he thesis 

The thesis is organized in the following chapters 

Chapter II : describes a mathematical model for 
predicting interference magnitudes and it’s extension to a 
nonlinear integer program formulation of FPP . The FPP is also 
formulated as a graph coloring problem. The nonlinear integer 
program formulation for evaluating the magnitude of a special 
type of interference problem, the intermodulation interference, 
is presented. 

Chapter III : Solution methodologies are developed based 
on the graph coloring algorithms, of which the sequential 
algorithms are reviewed in detail. A graph partitioning 
procedure is examined next, in the context of following a 
di V ide-and-conquer strategy of coloring the large graphs 
encountered in the problem. A heuristic method of merging the 
colored partitions without violating the coloring constraints is 
presented. Finally, a dual backtrack algorithm for solving the 
nonlinear integer program formulation of the i ntermodulat i on 
interference evaluation problem is discussed. 

Chapter IV : The implementation of computer programs to 
solve the FPP, using the solution methodologies of Chapter III, 
IS described. A solved example using the implemented system is 
presented for a problem segment taken from the Indian Railways 
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Microwave network. 

Appendix 1 : The technical details of the operation of 
the network is explained in a simplified manner. The propagation 
formulae for interference calculations are also given. The 
mechanisms by which the interference problems arise are explained 
at the end. 

Appendix II : The frequency planning activities currently 
going on in the Indian Railways for their microwave communication 
network are described. 



CHAPTER II 


vlATHEMATICAL DESCRIPTION OF THE FREQUENCY PLANNING PROBLEM 


In this chapter, the mathematical model of interference is 
resented first, which also draws upon the real-life data 
:"ovided by ”The Approach Papers for Frequency Planning for 
idian Railways Microwave Network*' [153. Next. two alternative 
ithematical formulations are described for the Frequency 
.anning Problem (FPP), namely, the Nonlinear 0-1 Integer Program 
)del and the Graph Coloring Problem model. Each formulation is 
icompanied by a literature survey of its domain. Finally, an 
iportant type of co-site constraint - the intermodulation 
iterference - is discussed in mathematical terms, based on the 
terature published in this area. 

1 A Mat^hematical model for Interf erence 

Consider a communication network as described below 

i = an index which uniquely identifies a link of the 

network, where 1 < t < n; 

F^ = frequency assigned to link i, 

where F^ takes discrete values from a given set F 


of f requenc i es ; 



= upper limit for interference level in link i, in the 
same units in vohich interference is measured; 

CA^ = set of neighbouring links of link i within the limits 
of co-ordination area of link i, 

where co-ordination area of a link is the area around 
it such that other links within this area can 
potentially interfere with its operation; 

f^ = a mathematical function which expresses the 
interference level in link i due to interference from 
various links i n CA . (in the same units in which 
interference is measured); 

The constraint on interference level in each link of the 
network is expressed as 

f . (F . , F . . . . F^) < U . V t ; 

i ^ j h h 

where F . . . . F, e CA . (2.1) 

j h L 


In general, the expression for f^ is not available in an explicit 
mathematical form. As per ^^The Approach Papers for Frequency 
Planning for Indian Railways Microwave Network” [1?], the 
expression for interference between two links i and j* , 
contains two distinct and additive terms, a . . and IRF. . , where 

3- . . is the RF -de coup 1 i nq between the two links and IRF. . is the 
Interference Reduction Factor which depends upon the frequencies 


and F^. assigned to the links i and j respectively. 
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2*1#1 RF-decoupling 

RF-decoup 1 i ng between two links indicates the extent to 
which operation of one link is free from interference due to 
operation of the other link. Its magnitude a^^ can be expressed 
as 


a 

RF* 


where 


0 


D 


a 

T 


a 

R 


Abi 


Aa 


6 


+ a 4- a 4- 

T R 


AG 4- Aa 4- S 


20 log 


D . 

■D — 

U)S 


difference be tween 


propagation loss in interfering and wanted signal 
paths ; 

angle attenuation at interfering transmitter; 
angle attenuation at interfered receiver; 
difference between antenna gains at wanted and 
interfering si gna 1 t ransmi 1 1 ers ; 

difference between attenuations at wanted and 
interfering signal transmitters; 

additional obstructions in interfering signal 
path . 


Higher the value of a is, for two links, larger is the extent 

RF 

to which the two links are free from mutual interference. A 

couple of points need to be noted about 


(i) It IS dependent upon the terrestrial layout of links 
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and sites of the network, as also upon the hardware 
installed at the various sites. 

<ii) It is independent of the frequency assigned to the 
1 i nks . 


2.1.2 InterferGnc© Reduction Factor 

The RF-decoupl i ng a , obtained for two links i and 
j, assumes that links i and J have the same frequency. This is 
the worst possible case of interference between links i and j. 
By assigning different frequencies and F^. . to links i and j 
respectively, the RF-decoup 1 i ng is further increased by the 
amount IRF. . and therefore, interference is reduced. That is why 
the term IRF is called Interference Reduction Factor. 


There is no single or simple expression for calculating 
IRF^^. for any two links i and j. Curves (graphs) are available 
which depict IRF. . as a nonlinear function of AF , ., where 




In contrast with a , it is to be noted that IRF. . 
depends upon the frequencies allocated to links i and j. Even 
the graphical curve for IRF^^. is not unique, the curve to be used 
depends upon the fixed characteristics of the radio relay systems 


used by the links i and j (FM 120, FM 500, FM 960 and ORS 34). 
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Since there is no easy expression for IRF. . 
calculations, following schemes may be considered 

(i> TaJbxilatton : Tables will be created for showing the 

variation of IRF. . with AF . .. Based on the fixed 
characteristics of the radio relay systems used by 
links t and j , a particular table may be selected and 
using AF . . = iF. - F .1 , the value of IRF. . may be read 
out . 

(ii> Piece-wise Linear Approximation of IRF^j cxirves : The 

piece-wise linear functions developed as an 
approximation to the original IRF curves may be used in 
their mathematical form to calculate IRF at any given 
AF. 

8.1.3 The Basic Interference Equation 

So the basic mathematical model of interference (2.1) 
now gets modified as in Equation (2.2) which follows. 
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IRF . . = f . . < If. - F .1 > , and 

^ J' 

a . . is independent of F. and F .. 

RF tj ^ t J 


2.2 A Nonlinear Integer Program formulation of FPP 


The basic nature of the Frequency Planning Problem 
(FPP) is allocation from a set F of frequencies, exactly one 
frequency k to each link i in the communication network. To 
formulate this problem , the decision variables chosen are 
assignment variables of 0-1 integer type, i.e. , 

Decision Varicible : 



1, if frequency k is assigned to 
link t ; 

0 , else. 


Constraints : 

Cl E X.. =1, Vi. where 1 < i < n. (2.3) 

A e F 

The above constraint implies that each link i is assigned exactly 
one frequency k from the set F. 

C2 D(i) > d^ , V i, where 1 < i < n; (2.4) 


D(i) = an expression for RF-decoup 1 i ng of link i from the rest of 
the links in the network as a result of a particular 
frequency allocation; and 

= minimum level of RF-decoup 1 i ng required for link t. 
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2.2.1 Expression for DCi) 

Given a table of IRF values which can be referred by 

the values IRF. . for any two links i and j having frequencies F. 

^ J? ^ 

and F^. respectively, or alternatively, given a function IRF^^. to 

calculate the IRF values, where 

IRF. . = f . . ( If. - F .1 >, 

' J- J ' 

D(i) can be evaluated in the manner described below. 


Assume frequency k is assigned to link i, and frequency 
L is assigned to link j, i.e. . 


Then , 


X.. = 1. 

zh 


and 


<F 


h'> 



1 . 

I ) 


E 

J € CA . 


E 

l € F 


IRF 


J i 


= 0(1 > , 


i f 


However, k is unknown apriori ; therefore 


ih 


1 . 


(F. = k) 


<2 . 5 > 


D(l ) 

E 

k € F 

< E 

J € CA^ 

E 

l € F 

IRF . . 

-it 

X 

y-S 

r.J 

X 

(2.6) 

This expression 

for D < i ) 

is to 

be read 

i n 

con J unct i on 

with 

constraint Cl so 

that summation over k (or 

1 > 

be 1 ong i ng 

to F , 


makes sense . 


Some points are to be noted in the expression for 0(i) 


as derived above 
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(i> It was assumed, by summing over links j in , that 

there is negligible interference fr om all links wh i c h 
are not in CA . . 

T. 

(ii) By summing over ail links j in CA^ , it was also assumed 
that RF-decoupi ing of link i from rest of the links in 
the network is the sum of RF-decoupl ing values of link 
i from each link j in CA.. 

There is another possibility of formulating the 
constraint C2 for every pair of links (i,j) where j is 
in CA^ , by considering individually the RF-decoupl i ng 
of link i with each link j in CA^ , instead of 
performing a sum over them as has been done above. 


(iii) The expression for D<i) is nonlinear as it contains 

product terms like <X.,.X.,). 

jl 

2.2.2 The Nonlinear Integer Program for FPP 


The nonlinear integer program formulation in it's final 
form is as given below. 

Cl E X., =1. Vi. where 1< i < n. (Z.7) 

k eF 


C2 


T. T 

M € F J- € CA. 

i. 


E 

l € F 


IRF . . X 


ih 



> d . . Vi. 


( 2 . 8 ) 
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where 1 < i < n, and 

X . , = 0 or 1 . 

\h 

It is noted that there is -no objective function in the 
above formulation. The objective function is trivial in this 
case because any feasible solution to the above problem is 
sufficient to provide a frequency plan for the network. Broadly 
speaking, the objective is to conserve the spectrum resources 
t.e. to use minimum number of frequencies in preparing a plan. 
The intermodulation phenomenon (Refer Appendix I for details) has 
not been considered in the above problem formulation, so the 
frequency plan obtained by solving the above problem is subject 
to verification for the presence or absence of i ntermodu lat i on 
interference. The problem of verification for the presence or 
absence of i ntermodu lat i on interference at any site in the 
network is formulated in a later section of this chapter. 

It is also noted that the problem size is very large 
for a communication network akin to the Indian Railways Microwave 
network. For allocating about 60 channels. available from 4 
re lay systems in use, to approximately 500 active links there 
will be 60 X 500 = 18000 decision variables of 0-1 integer 

type. The number of nonlinear constraints will be equal to the 
number of links i.e. about 500. Linearisation of the quadratic 
terms in equation (2.8), as a part of the solution strategy for 
solving the nonlinear program, will further multiply the number 
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of variables as well as the number of constraints. 

2.2.3 Literature Survey 

There are many approaches to solve a nonlinear integer 

program. Hansen C143 has reviewed the various methods of 

nonlinear 0-1 integer programming based on an extensive survey of 

existing literature. Three main approaches have been : 

enumerative methods, boolean algebraic methods and reduction to 

equivalent problems. Enumerative methods are basically branch 

and bound schemes of solving the 0-1 nonlinear program directly. 

Reduction to equivalent problems involve converting the problem 

to a polynomial 0-1 program for which many algorithms have been 

devised. One of the most important techniques is linearization 

in which the polynomial 0-1 program is reduced to a linear 0-1 

program by replacing all product terms FI x . by a new variable 

t € J 

z. Two new constraints get added : 

x^ + z+ jJj-l>0, and 

t € J 

- r - I j| z i 0. 

i ^ 3 

Even small nonlinear 0-1 programs can get enlarged considerably 
by this method, Glover and Woolsey [15] have proposed several 
ways of economizing on the number of new constraints introduced 
in the process of linearization. A more recent survey of the 
various general algorithms developed for a nonlinear integer 
program appears in Cooper [16]. 
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2. 3 A Graph Coloring Problem Formulation of FPP 

A Graph Coloring Problem formulation of the Frequency 
Planning Problem can be obtained by viewing the commun i cat i on 
network as a graph G<N,A>, where 

N = set of sites in the conrimun i cat i on network; and 

A = set of links where each link joins two sites in N 

un i que 1 y . 

Based on the geographical location of the links (and sites) and 
the communications hardware assemblage at the sites ulz. antenna, 
transmitter/receiver. bandpass filters etc., link-to-link 
interference potential (IP) is calculated for every pair of links 
within co-ordination distance of each other assuming that all 
links are operating at the same frequency. These values are 
stored in a ma t r i x IP, 

IP = (IP. ,) , where ( t , j ) e A , and j' € CA . or i e CA , . 

Stipulating a tolerance level of this interference potential for 
each link i (based on the number of links in CA^ > , a secondary 
graph GC is developed in the following manner 

(i) Each link i in the graph G becomes a vertex i in GC, 

(li) Two vertices i and j in GC are joined by an edge if the 
two links i and j they represent have a mutual 
interference potential greater than the stipulated 
tolerance level for either link i or link j. 



2 ? 


The 

secondary graph 

GC 

i nd i cates 

al 1 

the critical 

interf erence 

linkages between 

the 

links of 

the 

commun i cat i on 


network. The criticality here is defined with respect tp the 
stipulated level of tolerance of the mutual interference 
potential between two links. 

The objective of frequency planning is to allocate a 
discrete frequency to each vertex of GC which is a link of the 
communication network. If two vertices which are joined by an 
edge in GC are assigned the same frequency then the interference 
between the two links they represent will obviously be 
unacceptably high because it exceeds the stipulated tolerance 
level of interference (This is nothing but the Co-channel 
constraint described earlier in Chapter I and explained in in 
greater detail in Appendix I>. However if they are assigned 
different frequencies then the interference decreases sharply and 
there is no interference problem of the original order of 
magnitude. Since the frequencies which are allocated are drawn 
from a limited set F, it makes sense to minimise the actual 
number of frequencies used in the allocation process. It might 
also be desirable to conserve the scarcely available frequencies 
for fulfilling the future requirements of the communication 
network . 


By treating the discrete frequencies as colors. it is 


now possible to visualise the FPP as a Graph Coloring Problem on 
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the graph GC . In a Graph Coloring Problem, it is sought to color 
all the vertices of a graph using minimum number of polors such 
that no two adjacent vertices (vertices which are joined by an 
edge) have the same color. 

/ 

2.3.1 Additional constraints on Graph Coloring based on 

Adjacent channel constraints in Frequency Planning 

Sometimes, even if two neighbouring vertices in GC are 
assigned adjacent frequencies (adjacent in the frequency 
spectrum) there might be interference problems (This is how 
Ad jacent -channe 1 constraints arise). So the frequency allocation 
constraint is modified in the sense that two neighbouring 
vertices in GC cannot be assigned frequencies that are the same 
or adjacent to each other t.©., they must have a separation of at 
least one frequency channel. For the Graph Coloring Problem. it 
IS equivalent to saying that all the vertices of the graph GC are 
to be colored using minimum number of colors such that no two 
neighbouring vertices are assigned the same or adjacent colors. 
Two colors are adjacent if the frequencies they represent are 
adjacent in the spectrum. 

2.3.2 Some charactorislics of th© f orimilalion 

Following points may be noted in the formulation 


described above 
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(i> The above process of tightening the frequency 

allocation constraint may be extended further. For 
example, the frequency allocation can be constrained 
such that two adjacent vertices in GC cannot be 

assigned frequencies that are separated by less than 2 
< or 5 or 4 or . . . ) frequency channels. 

(ii> Once again, intermodulation interference is not 

considered in the above formulation. It has to be 
checked separately after the frequency allocation is 
specified completely. As stated earlier, the 

formulation of the problem of evaluating 
i nte rmodu lat 1 on interference is taken up under Section 
4 of this chapter. 

2.3.3 Literature Survey of the Graph Coloring Problem 

The Graph Coloring Problem has been widely researched, 
mainly on account of the famous Four Color Problem (FCP). The 
Four Color Problem is the problem of coloring a planar graph 
using at most 4 colors. This problem was posed more than a 
hundred years ago by cartographers (map-makers) who found, by 
experience that they could color any map, in which any two 
countries sharing a common border are colored differently, with a 
maximum of 4 colors (A geographical map can be represented by a 
planar graph). But no proof of this problem was forthcoming 
until a decade back when a computer was used to analyse around 
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thousand special cases of the FCP to shoiA» conclusively that the 
FCP had indeed a proof and a solution. This solution is due to 
K. Appel, IaI, Haken, and J.A. Koch [17]. It may be noted that the 
graph GC , as derived from the communication network, is rarely 
planar in nature. 

In light of the fact that the FCP is a special case of 
the Graph Coloring Problem, and has itself defied a solution for 
close to 1^0 years - being solved ultimately only with the aid of 
a computer, it would seem likely that there is something 
inherently difficult about the Grapfi Coloring Problem. Indeed 
with the development of the theory of computational complexity of 
combinatorial problems, mainly by Garey and Johnson [183, it has 
been found that the Graph Coloring Problem, also known as the 
Chromatic Number Problem, is NP-Complete ([193, [20]). 
Therefore, it seems unlikely that any polynomial time-bounded 
algorithm exists for solving this problem optimally. 

Graph coloring has considerable applications in the 
areas of combinatorial optimization. Problems involving conflict 
resolution or optimal partitioning of mutually exclusive events 
can be modelled as graph coloring problems. The examination 
s c e d u 1 1 n g problem, in which ex am i nations are to be scheduled in 
the smallest number of time periods such that no individual is 
required to participate in two examinations simultaneously, is 
modelled as a graph coloring problem with examinations as 


vertices of a graph and time periods as colors. 


Two vertices are 
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joined by an edge if both the corresponding examinations require 
the par t i c i pat i on of an individual. So the problem of coloring 
the graph using minimum number of colors is equivalent to using 
minimum number of time periods to schedule the exami na t i ons . 

Storage of chemicals on the m i n i nnjm number uf shelves 
such that, no two mutually dangerous (reactive) chemicals are 
stored on the same shelf, can also he modelled as a graph 
coloring problem. Machine and job scf'ieduling and loading 
problems have also been defined as graph coloring problems. 

The constraints in a problem, which can be modelled as 
a graph coloring problem, are usually expressible in the form of 
pairs of ( 1 n ) compa t i b 1 e objects. Such ( i n > compat i b i 1 i t i es are 
represented by the edges of the graph where the objects become 
the vertices. In practice, the constraints imposed are more 
complex than i ncompat i b i 1 i t y constraints ( Ad jacent -channe 1 
constraint is one example) and hence the models developed are no 
longer simple coloring problems. 

The methodologies known for the solution of the Graph 
Coloring Problem are reviewed from published literature in the 
first part of Chapter III. 
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2# 4 A Ma-thematical Etescription of InternK>dulation Interference 

An important type of co-site constraint is the 
int&rrr^J^d^ilation inteTf&T&rvce which may be expressed 
mathematically unlike some other constraints on frequency 
planning. The intermodulation phenomenon is the creation of 
secondary frequencies in the transmitter/receiver circuitry as a 
result of interaction of multiple active frequencies at a site. 
These secondary frequencies can interfere with other signals and 
distort the quality of radio communication. Thus in a large and 
dense commun i cat i on network, it is very crucial to minimize this 
type of interference even though it may be of secondary nature, 

2# 4# 1 Mathematical condition for the occurence of 

Inter modulation Interference 

Consider a site which is communicating directly with n 

other sites using n links. Let be the frequency assigned to 

link I, where 1 < i < n. The intermodulation interference 

evaluation problem is to determine if the set of frequencies , 

F., .... F operating simultaneously at the site is free of the 

L n 

intermodulation phenomenon upto a pre-spec i f i ed order Q or less 

(Order Q is defined a little later). It turns out that the 

transmission frequencies create secondary frequencies, which are 

integer linear combinations of F,, F..., F . If one of these 

1 Z n 

secondary frequencies is the same as a receiving frequency, which 

could be one of the frequencies F,, F^, .... F . then 

12 n 
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i nte rmodu 1 at 

ion interference 

occurs. It should 

be 

no t ed 

that 

each of the 

frequencies F^ , 

F , .... F can 

2 n 

be 

used 

as 

a 

t ransmi ss i on 

frequency only 

or a receiving frequency 

only 

at 

a 


given time. Mathematically, intermodulation interference 

therefore occurs whenever 


n 


E F. X, = 

t = 1 

0. 

for some 

integers X^ where 

1 

< t < n , 


and at least one 

of the 

X . takes 

the value ± 1 . 




The sum 

Q = 

n 

E hj 

is def i ned as 

the 

order 

of 



i=l 





interference. It 

turns 

out that 

the interference 

i s 

strong 

if 

the order Q is 

smal 1 . 

Hence 

finding the lowest 

order 

of 

i ntermodu lat i on 

is useful in 

evaluating a given 

frequency 


assignment, in this case, the set of frequencies F,. F., .... F 

1 Z n 

assigned to the links l..n incident at the site under 
cons 1 derat i on . 

2.4.2 A nonlinear integer program formulation for evaluation 

of Inter modulation Interference 

The problem of evaluating whether the frequency 
assignment to the links incident at a site is free of 
intermodulation interference has been formulated by Mathur et al 
[8] and can be stated as 

n 

Mm. = E <2. 9) 
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n 

s.t. r F, X. = 0. (2.10) 

i = l 

where 1 < t < n, and X are integers, and 
3 at least one i s.t. | 1 = 1- <2. 11) 

opt 

If Q„ > Q + 1 then the set of frequencies F, . F^ F can 

u 1 2 n 

be declared as i nt ermodu lat i on-free . 

In the real world, the receiving frequency is 
interfered with if the secondary frequency due to intermodulation 
IS close, but not necessarily equal, to the receiving frequency. 
As a result the constraint <2.10> is modified so that 

n 

-GB < X! F^ X^ < GB, (2.12) 

11=1 

where GB i s a known parameter called Guard Band. 

It may be noted that the above problem has to be solved 

for each site in the network. Moreover it can be solved only 

when the frequency assignment is known aprior-i for the 
communication network. A solution to this problem does not 
generate a frequency plan, rather it evaluates the given 

frequency plan with respect to the phenomenon of intermodulation 

interference at the various sites in the network. However this 
will aid the process of preparing a feasible frequency plan. 
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2.4.3 Lit.erat.ure Survey of solution methodologies 

From the problem formulation, it is clear that it can 
be reduced to a nonlinear integer program, the modulus operator 
in the objective function being the source of nonlinearity. 
While any of the general algorithms for a nonlinear integer 
program may be used for solving the above problem, efficient 
solution methods must be developed, as this problem has to be 
solved many number of times for evaluating a given frequency 
plan . 


Two backtrack algorithms proposed by Morito et al C213 
have been examined by Morito et al C22], one solves the primal 
problem while the second solves the dual. The dual backtrack 
algorithm appears to be computationally favourable and is 
discussed in greater detail in Chapter III. Nishimura et al [253 
also investigate two branch-and-search techniques. again for 
primal and dual versions of the problem, for generating subsets 
of frequencies which are free from intermodulation. 

2.5 Final Choice of the Problem Formulation 

The Graph Coloring Problem formulation is taken up for 
developing solution methodologies for the FPP , in preference to 
the Nonlinear 0-1 Program approach. The latter formulation is 
very bulky in size and is not handled easily from computational 
point of view. The graph coloring approach has a compact 
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structure in the form of graph and it has a better defined 
objective of minimizing the number of colors used i.e. ocnserving 
the frequency resources. Moreover, a large variety of graph 
coloring algorithms are available in literature for ready use. 



CHAPTER III 


SOLUTION METHODOLOGY FOR THE FREQUENCY PLANNING PROBLEM 


In this chapter, the solution methodologies for solving 
the FPP are developed in accordance with the Graph Coloring 
Problem formulation, described in the previous chapter. 

Initially the various ^rctph. coloring ro-atin&s are explained 
followed by a description of the modifications necessitated in 
order to accomodate the Adjacent-channel constraint in frequency 
planning. As the graph GC, derived from the communication 
network, is large-sized and the interference effects are of local 
nature, a diivide-and-conqxier strategy is devised for partitioning 
thje graph based on it's edge connectivity. After recursively 
coloring the partitions, a procedure is developed for merging 
them back to the original graph GC, reconciling the vertex 
colorings wherever any conflict arises. All the procedures 
mentioned above are presented along with the basic theoretical 
background. Additionally, reference is invited to the literature 
published in the associated domains. The solution methodology 
for carrying out the intermoduLation interference evaluation of 
the frequency plan generated by the coloring of the graph GC is 
also discussed based on the nonlinear integer program model given 
under section 4 of Chapter II. 
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3.1 Graph Coloring Algorithms 

There have been two major approaches to vertex coloring 
of a graph - irude pendent set approach and seqxLential approach. 
As is to be expected, each approach has it's own advantages and 
drawbacks. Of course, there are also hybrid approaches which try 
to seek the benefits of both the approaches by combining the 
strategy of sequential coloring with independent set approach. 
Both the approaches are discussed below, under separate sections. 

3.1.1 Independent Set Approach 

Some definitions borrowed from Syslo et aL C241 
are necessary in order to explain the independent set approach. 

[Definition! A subset of vertices of a graph in which no 

two vertices are joined by an edge is called an independent set .M 


[Definition! A subset W of vertices is a 

independent set (MIS) if there is no independent set in 

is contained properly. 

[Definition! A maxirn-um. independent set (MmIS), on the 

other hand, is a maximum cardinality independent set. ■ 

A coloring of the vertices of a graph G using k 
colors is equivalent to a partition of the set of vertices of G 


maximal 
which W 



into k independent sets V,, , .... V, . such that 

1 Z M 

n ^ ^ ^ ^ J> 3Lnd i,J = 1 k, and 

h 

U V. = V. 
i=l ^ 

Such a partition is called a k-coloring partition of V. The 

independent set approach is to discover such a partition by first 
coloring vertices with color 1 i.&. finding set , then set . 
and so on upto set V^. 

To understand the basic relationship, following which 
various algorithms, both exact and approximate. have been 
developed for solving the graph coloring problem, Theorem (5.1) 
is quoted below from Syslo et al C24]. 

[Theorem 3.11 For every /^-coloring partition 

# ' I 

of a graph G there exists a l-coloring partition ^ 2 » 

t 

> of G such that I < h and at least one of the sets in the 
latter partition is a MIS. 

Proof If Vj is not a MIS of G then it can be augmented 

t 

by the vertices of the other subsets to become such a set . 

« f 

Now. every set ~ ~ ^ ^ 

i ndependent , however , some of them may be empty. Dropping empty 

1 f 

sets = <p , if any, form a^ l-coloring partition of G for i = 

f 

1. 2, .... k where I < K, and is a MIS. ■ 

It follows from Theorem <3.1) that in every partition 
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of G due to an optimal coloring, all sets are MIS. Let W be one 
such MIS of G. Then all vertices in W can be colored by the same 
color. Let be the subgraph induced by the set of vertices 
V - W. Hence 

;t<G) = Min + 1 (3.1) 
W c V 

where : 5 :<G) denotes the number of colors used in an optimal 
coloring of G, also called the chromatic number of G and it is 
assumed that ;^<0) = 0. 

Equation (5.1) is the basis for the various algorithms 
which follow the independent set approach. 

3.1. 1.1 An exact method for optimal vertex coloring of a graph 

An exact procedure for coloring a graph using the 
independent set approach is a backtracking method, as most of the 
other exact methods are. If a graph G is >?-chromat i c , then it 
can be colored with k colors; coloring first with color 1 a MIS 
W, of G, next coloring with another color 2 a MIS W„ of G|,, ,, , 

1 2 'v-Wj 

and so on until all the vertices of G are colored. However, this 
method does not specify which MIS should be colored at each step; 
therfore, every MIS of a current subgraph must be considered. 
That is precisely the reason for backtracking in such a method. 
An algorithm due to Christofides C253, illustrates the above 


procedure . 
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3. 1.1.2 Approximate algorithms 

The idea of generating a set of vertices of one color 
at a time gave rise to several approximate algorithms. The 
maximum independent set algorithm (MmIS algorithm) generates a 
sequence of maximum cardinality independent subsets in subgraphs 
of uncolored vertices and colors them one by one. 

The subproblem of finding a MmIS in a graph has no 
po 1 ynomia 1 -t ime algorithm, so an efficient algorithm is used to 
generate an approximately MmIS. This becomes the genesis of the 
approximately maximum independent set (AMmIS) algorithm. 

Both the versions of the approximate coloring algorithm 
can be found in Sysio et al C24]. Interestingly enough, while 
the above approximate algorithms, based on the independent set 
approach, exhibfit better worst-case behaviour than the algorithms 
based on the sequential approach described below, in their 
average behaviour they produce substantially inferior colorings. 

3.1.2 Sequential Approach 

Let V,, V- V be an ordering of the vertices 

12 n 

of a graph G. In a sequential method for coloring the vertices 

of G, vertex v^ is added to the subgraph induced by already 

colored vertices v, , v., v. . and a new coloring of v,, v„, 

IZ i-l IZ 

.... determined . This step is repeated for i = 1, 2, 
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.... n, where for i = 1 the subgraph is empty. At each step an 
attempt is made to use relatively small number of colors. 

The following sequential algorithms are reviewed 
as in a paper by Matula &t ctl C263. 

3. 1 . 2. 1 Sinq^ld Sequent.ial Algorithms 

« 

A simple sequential algorithm goes as follows : - 

ALGORITHM S IMPLE_SEQUENT I AL ; 

: COLOR(Vj) := 1; 

: FOR i := 2, 5 n DO 

COLOR(w.) := MIN C (k : k > 1> AND COLOR < .) k 

^ J 

for every v^. (1 < J < i') 
adjacent to v^ > 

It is noted that in this simple algorithm the vertices 
colored prior to coloring v^ retain their colors. Each vertex v^ 
can be colored by color i, therefore COLOR<v^> < i. On the other 
hand, at least one of the first degCv.) + 1 colors can be 

assigned to v.. Hence, 

CDLORlv^) < Min Ct , deg<v^) + 1> (5.2) 

for every i - 1, 2, .... n, and thus 


STEP 1 
STEP 2 


END. 
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X <G) < Max <Min £i . deg(v.> + 1>>. (5.5) 

S 1 < £ < n ^ 

A1 ternat i ve 1 y . 

u»<G; V,, Vj, .... V ) = Max (Min ££ , deg(v.) + 1>) (5.4) 

1 < L < n 

is the upper bound on the number of colors used by the Sequential 

algorithm to color the vertices v, , v,. v of graph G. This 

1 z n 

holds true irrespective of the, ordering of the vertices of the 
graph G. 


There are various orderings of the vertices of G which 

have been considered for applying the Sequential algorithm. An 

obvious beginning was made with the LF ordering in which vertices 

with the largest degrees were ranked first in the order i.e. 

deg(v,) > deg(v„) > ... > deg(v ). The sequential algorithm run 
14 n 

on such an ordering is called the largest first secfuentiaL 
algorithm (LFS algorithm). An LF ordering is superior to the 
other orderings when they are used in the sequential algorithm. 
This fact is stated as Theorem (5.2) below, again quoted from 
Sys lo et ctl [24] . 


C Theorem 3. 21 


. Let u, , u^, . . . , u be an -LF ordering of the 

14 n 


vertices of a graph G and u_(G; v, , v^ v ) be the upper 

b 1 4 n 

bound on the number of colors used by the Sequential algorithm 
for coloring G with any arbitrary ordering of it's vertices. Let 


u, r- = u,.^(G; u, , u„ u ) be the upper bound on the number of 

Lr b 1 4 n 
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colors used by the Sequential algorithm for coloring G with an LF 

ordering u^ of it's vertices. Then 

u = Min (G; n, . n„. .... n > 

LF 1 2 n 

where the minimum is taken over all orderings n, , n*. .... n of 

12 n 

the vertices of G. ■ 

(The proof, which is omitted here, is outlined in Syslo et al 
C24] . ) 

Another ordering which is similar to LF ordering is the 
SL <.sma.L lest-last } ordering. Such an ordering is constructed for 
the vertices of graph G using the following procedure described 
by Matula C26] 

<i> v^ is a minimum degree vertex of G. 

(ii) For i = n-1 , n-2, .... 2, 1, 

is a minimum degree vertex in the subgraph of G 

induced by V - Cv , v ,, .... v.,,}. 

n n-1 L+1 

The ordering is called smal 1 est - last because 

deg.(v.> = Min Cdeg.<v.>} (5.5) 

1 < j < i 

where deg^(v^) is the degree of v^ in the subgraph induced by v^^ , 

v^ . The sequential coloring algorithm applied to a SL 

ordering is called the SLS algorithm. 
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A careful examination of the Sequential Coloring 
algorithm reveals that a vertex v^ can be colored by one of the 
first i colors or by one of the first deg.(v.> + 1 colors, 

whichever is minimum. Hence the upper bound v^ . v^ 

v^) on the number of colors used by the Sequential algorithm can 
be improved to 

u«'(G; V,, V-. .... V > = Max Min Ct . deg.<v.> + !> 

1 < t<n 

(5.6) 

Since deg,(v.) < deg<v.> for any ordering v, , v-, .... v of the 

t t t 1 2 ^ n 

vertices, therefore 

u„' (G; V . V- V ) < Uo(G; V,. v, V >. (5.7) 

bl2 n S12 n 

Let u». = u_'<G: t, , t,, t ) be the upper bound on the 

oL b i c n 

number of colors used by the Sequential algorithm for coloring G 

with an SL-ordering of it’s vertices t,, t „ , .... t . By the 

12 n 

nature of construction of SL-ordering of the vertices. 


u„, = Min u,.’<G; n, . n, n >. (5.8) 

SL 5 12 n 

where the minimum is taken over all the orderings n, , FI. FI 

12 n 

of the vertices of G. 

By (5.7), u„'(G; u.. u, u_) < u_(G; u, , u„ u ), 

S 12 nS12 n 


(5.9) 


where u, , u_ u is a LF-ordering of the vertices of G. 

12 n 


u„’(G; t , t t ) < u ' (G; u , u,, ..., u ) 

S 12 nb 12 n 


By (5.8). 
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<3. 10> 


Therefore , 


^SL ^1’ ^2’ '^1' ^2' ■'■’ *^0^ ~ '^LF’ 


i .e. , 


^SL - ^LF 


(3.11) 


which suggests that the SLS algorithm is superior to the LFS 
algor i thm. 


There are several other heuristic algorithms for graph 
coloring which are sequential in nature. An important one 
amongst these is the SLF (satxiration-largest~/irst ) algorithm, 
also known as the DSATUR algorithm, due to Br61az C27]. 

IDefinition] A satviration degree of an uncoiored vertex v 

in a partially colored graph G is defined as the number of 
different colors to which v is adjacent. ■ 

The SLF algorithm proceeds as follows 

(i> Color a vertex of minimum degree in graph G with color 
1 . 

(ii) For i = 2, 3, .... n do 

As V. take an uncolored vertex of maximum saturation 

3 , 

degree breaking ties, if any, by choosing the vertex of 
greater degree and color v^ with the least possible 
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color. 


3.1. 2.3 Sequential Algorithms with Interchange CSID 

The simple sequential algorithm can be improved by 
re-arranging, at each step in the coloring sequence. the colors 
of the already colored vertices to avoid, if possible, using an 
additional color for the next vertex to be colored. The 
interchange procedure in Sequential algorithms was first 
described by Matula et al C26]. 

Assume that , V 2 been colored by k 

colors and the simple sequential algorithm colors v^ with color 
k+1 . This implies that v^ has neighbours amongst the vertices 
Vj , v^. .... colored with 1, 2. .... h. If the subgraph 

formed by these neighbours of v^ contains a complete subgraph on 
k vertices then the new color is necessary for v^ . 

Otherwise, it may be possible to interchange the colors of some 
neighbours of v^ , preserve the As-coloring of the subgraph 
generated by v^ , v^, .... free one of the first k colors 

for V . . 

I. 

Let G be the subgraph induced by the vertices colored 

fXf 

with color p or <?, 1 < p, g < M. If G has two vertices u and v 

pcf 

of different colors, both adjacent to v^ and both belonging to 

the same connected component of G , then no interchange of p and 

P9 

g can free one of these two colors for v.. Therefore, if this 
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holds for every p and q, must be colored with M+1 . On the 
other hand, if there exist p and q such that in every connected 
component of vertices adjacent to v^ are of at most one 

color, then the <p,<?>- interchange procedure, described below, can 
be applied fruitfully. 

In the <p ,c?) - i nterchange procedure let U be the set of 

P 

vertices of those connected components of G which contain a 

P<? 

neighbour v of vertex v^ , such that COLORCv) = p. Now, for every 

u € U , if COLORiu} = p then set COLORiu) = q, and if COLOR(u> = 
P 

q then set COLORiu) = p. 

It is easily seen that while the coloring of v^. v^. 

..., V. , remains feasible, color p has become free for v.. The 
interchange procedure, as described above, is an extension of the 
original version due to Matula et al CZ63 and it initially 

appeared in a paper by Johnson C28]. 

Obviously the SI algorithms can be applied to any 

ordering of the vertices of a graph, for example one can have 

LFSI and SLS I algorithms. 

Despite the claim that the SI algorithm is an 

improvement over the Simple Sequential algorithm, it is possible 
that the SI algorithm uses more colors than the Simple Sequential 
algorithm for some random graphs and orderings. This is 
particularly true if the interchange of colors is carried out 
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early in the sequential coloring procedure. In general, all the 
approximate coloring algorithms following the sequential approach 
are capable of producing very bad colorings, in terms of number 
of colors used, for certain graphs. On the average, however, 
they seem to fare reasonably well, both in terms of running time 
as well as number of colors used. 

3.1. 2.3 Sequent-ial Algorithm with Backtracking C SB ) 

The idea of sequential coloring can be extended to 
give rise to a backtracking method which will color every graph 
optimally i.e. with the minimum number of colors. The 
backtracking sequential algorithm was first proposed by Brown 
C29] . 


Consider a graph G with an arbitrary ordering of it's 

vertices v, , v., v . Let v, , v. v. , be vertices 

I c n 1 Z t-l 

which are already colored. For the vertex v^ , find a set of 

feasible colors and assign Min CU.} to vertex v.. Let 1. , 

denote the maximum color number used to color the vertices v^^ , 

v^ Every color J that belongs to must satisfy 

each of the following conditions 

(i) j ^ ^i-1 ^ that no redundant colorings 

are generated at any stage. 

Min Ci, deg(v^> + 1>, 


( 1 i ) J 


as per < 5 . 1 ) 
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<iii> J is not a color of any vertex (1 < h < t-1) which 
is a neighbour of . 

(iv) If a complete qi-coloring of G is known, then j < q - 1. 

The coloring procedure backtracks from vertex to 
vertex whenever ^ at any stage. This is done by 
deleting from U. , the current color of v . , and restarting the 
coloring procedure from the vertex v^ ^ using the new The 
procedure exits if a backtrack is required from vertex v^ 
implying that no coloring better than the already known 
(^-coloring can be found. The procedure also exits when is the 
last vertex colored and this implies that a l-colonng has been 
found, where I < q. Set q = I and try to get a fresh coloring by 
repeating the backtracking procedure afresh. 

The backtracking algorithm described above may be 
refined by using an intelligent ordering of the vertices. If 
instead of an arbitrary ordering, the first h vertices are so 
ordered that they form a clique (each of the vertices in a clique 


is adjacent 

to all the others ) 

, which 

i s 

the largest in 

G. 

f o i lowed by 

the renrtaining 

vertices 

V . in 

1 

such 

an ordering 

such 

that V . IS 

L 

adjacent to more of 

the 

vertices 

'^i- 

^t-l 

than any other vertex v . 

J 

for J 

> i 

, then 

the 

backtrack i ng 

need 

not go back 

farther than 

'^h+1 * 


In other 

words. it is 

not 

necessary 

to backtrack 

to 

''l 

for 

enumerat i ng all 

the 
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non-redundant colorings. So the algorithm becomes more 
efficient. However, finding the largest clique in a graph is 
also a NP-Complete problem like the Graph Coloring Problem; hence 
a greedy method of ordering the vertices due to Brown C293 is 
chosen. In this ordering, called the Greedy Largest First 
Ordering (GLF) . the first vertex is chosen to be of the 
largest degree in G. Subsequent vertices v^ are such that is 

adjacent to more of the vertices v^, ^'^7 other 

vertex v. forj' > i . 

Look-ahead procedures can also be incorporated into the 
backtracking algorithm to improve it's efficiency. These look a 
few steps ahead and check the effect of a current color 
assignment on a later stage and if any infeasibility is detected 
early on, the current color assignment can be dropped and the 
next color can be considered for assignment right away. Before 
applying such methods, the additional costs of computation for a 
look-ahead procedure must be balanced against the savings in 
enumeration it generates for the backtracking algorithm. The 
look-ahead procedures for sequential backtracking algorithm have 
been explored by Brown 1291 and Korman [30] . 


3.2 Partitioning the graph into easily colorable graphs 

The technique of solving a problem by breaking it into 
smaller, more easily solvable parts, solving these partial 
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problems, and then combining the partial solutions into a 
solution to the whole problem is called Di-oide-an-d-Conquer. 
Applied to the design of algorithms, the method involves 
partitioning a problem into smaller, homogenous pieces, to which 
the overall technique is then applied recursively. If the parts 
can be combined efficiently, the method may lead to an efficient 
algorithm. PartitioniTxg is a related technique, in which a 
problem is divided into pieces, solutions are obtained for each 
piece, and are then combined to form a solution for the whole 
problem. The rather thin difference between di vide-and-conquer 
and partitioning is that d i v i de-and-conquer typically involves 
recursion, while partitioning does not. However, this difference 
will be ignored in the following text. 

The di V ide-and-conquer technique is used in many graph 
algorithms. To apply this technique, the input graph is 
decomposed into a hierarchy of components. Then the problem is 
solved on each of the components at the bottom of the hierarchy 
and gradually the solutions on larger and larger components are 
pieced together, until a solution is finally found for the entire 
graph . 


Usually the decomposition that supports 
d iv i de-and-conquer is based on finding separators of the smallest 
size. Separators are those entities of the graph whose 
elimination results in the decomposition of the graph. Different 
kinds of separators can be considered for decomposing a graph :- 
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(i) Vertex set - elimination of a set of vertices from the 

graph partitions it in two or more components. 

< F i gure 5 . 1 > 

(ii> Edge set - elimination of a set of edges results in a 

disconnection of the graph into two or more 

components. This is equivalent to generating a 

s-t cut in the graph, where s and t are vertices 
in the two different components. (Figure 3.2) 

<iii)Clique - elimination of a clique results in the 
decomposition of the graph into two or more 

components. (Figure 5.5) 

The separator chosen for the partition of the input 

graph into easily colorable graphs is the Edge set. The basis 

for this choice lies in the fact that the interference phenomena 
is of a local nature and the edges, in the graph GC to be 
colored, are representative of the constraints on coloring 
(frequency assignment) due to the interference effect. So by 
partitioning the graph GC while choosing Edge sets of least 
cardinality, the induced subgraphs are inter-connected by minimal 
number of edges and hence have the maximum de-coupling possible. 
In such a situation. the effect of coloring (frequency 
assignment) within an induced subgraph will have minimal effect 
over the coloring (frequency plan) for the other induced 
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FIGURE 3.1 Vertex separator 
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FIGURE 3.2 Edge separators 




FIGURE 3.3 Clique separator 
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subgraphs. As a result the amount of computational work required 
to combine the individual colorings of_all the induced subgraphs 
and get the coloring for the entire input graph, is kept as low 
as pos s i bl e . 


The procedure for generating the Edge set is based on 
the Max Flow - Min Cut Theorem due to Ford and Fulkerson C31]. 
All the edges in the input graph GC are assumed to have unit 
capacity and so the capacity of the minimum cut is equal to the 
number of edges present in the minimum cut. The edges in the 
minimum cut comprise the Edge set. 

[Definition] EdLge conn&ctivity between two subgraphs 

and . induced by two disjoint subsets and of the vertices 
of a graph G, is the total number of edges in G joining a vertex 
in Gj to a vertex in G^ . ■ 

The above definition implies that, given two subgraphs 
Gj and G^ induced by two disjoint subsets of vertices such that 
the edge connectivity between them is K, the number of edges that 
have to be cut in order to disconnect any two vertices - one in 
G^ and the other in G^ “ i s K or less. Given a specified maximum 
value of the edge connectivity K between any two subgraphs 
induced by the partition, it is desired that the maximum flow (or 
minimum cut) value between any two vertices in the same induced 
subgraph exceeds the value K i.e>. to disconnect any two vertices 
in the same induced subgraph, more than K edges have to be cut. 
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The partitioning procedure relies upon the knowledge of maximum 
flow values for all pairs of vertices in the original graph, 
which in turn can be computed efficiently by constructing a 
Gomory-Hu Cut-Tree C52]. 

2* 1 t^erminal maximum flovfs and Graph Partitioning 

The Gomory-Hu cut-tree which is used to compute mu 1 1 i -t ermi nai 
maximum flows i.e. maximum flow values for more than one pair of 
vertices, has the following properties 

(i) It is flow-equivalent to the original graph from which 
it is derived. Two networks (graphs or trees) are 
called flow-equivalent if there is one-to-one relation 
between the vertices of the two networks and the 
maximum flow values between any two of the vertices is 
the same in both the networks. 

(ii) Each link of the tree represents a minimum cut of the 
original graph. This is also the reason why it was 
called a cut-tree by Gomory and Hu. 

(lii)The complete Gomory-Hu cut-tree can be constructed by 
solving 0<n) maximum flow problems where n is the 
number of vertices in the graph. (The process of 
construction of a complete Gomory-Hu cut-tree is 


described under Section 2.2) 




A given graph GC 



Equivalent tree 

network which is a 
Gomory-Hu cut-tree 


FIGURE 3.4 Flow equivalence of a given graph to it’s Gomory-Hu cut -tree 
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(Figure 5.4 illustrates the Gomory-Hu cut-tree derived from a 
given graph and the flow equivalence of the two.) 

The maximum flow value between any two vertices is equal to 

the capacity of the unique path between the two vertices in an 

equivalent tree network. Using a predecessor index 

representation of the equivalent tree network, the path between 

any two vertices can be determined quite efficiently and hence, 

the maximum flow between the two vertices can also be computed 

easily. Therefore, the maximum flow values for all pairs of 

2 

vertices (there will be 0(n > pairs) can be found out by solving 

2 

0(n) maximum flow problems, instead of 0(n ). by virtue of 
property (iii) stated above. Given the maximum flow values for 
all pairs of vertices, the graph may be partitioned along 
following lines : - 

ALGORITHM PARTITION 
STEP 1 : I := 1 

CONSIDER a vertex v from V /*V is the vertex 
set of the graph to be partitioned */ 
Pj := Cv) /* Pj is the I partition subset- 

of V 

V := V - Cv) 

READiK) /* K is the specified edge 

connectivity between any two partitions ■*•/ 


STEP 2 


IF \) <p THEN 



5 ? 


CONSIDER a vertex u from V 
J := 1 

STEP 2.1 : IF J < I THEN 

IF max imum-f low(u , Pj) > K for every 

vertex Pj e Pj THEN Pj : = Pj u Cul 
ELSE J := J + 1 
GOTO / 2.1 
ELSE I := I + 1 
Pj := Cul 
GOTO / 2 

STEP 5 ; FOR J ;= 1 TO I OUTPUTiPj} Pj is the 

partition subset of V, the vertex set of the 
graph */ 

3.2.2 Construction of a Gomory— Hu cut— tre© 

The partitioning algorithm described in Section 2.1 
utilized the maximum flow values between all pairs of vertices, 
which in turn could be calculated efficiently by using the 
complete Gomory-Hu cut-tree. The process of construction of a 
complete Gomory-Hu cut-tree (a few stages are shown in Figure 
3.5) involves the following steps :- 

STEP CiO Select two vertices - SOURCE and SINK - arbitrarily in 
the original graph and compute the maximum flow between 


them using the original graph. 
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,) First partition of the graph into 
bignodes X and X. Maximum 
flow between SOURCE and 

SINK is x^. 


if 



(b) Second partition of the graph by 
splitting bignode X into and . 



FIGURE 3.5 Stages in the construction of a Gomory- Hu cut- tree 
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STEP CiiJ> Use the resultant minimum cut to partition the vertices 
into two parts, SOURCE and SINK being in different 
partitions. Label each partition as a Bignode along 
with a unique index number attached to identify any 
b i gnode . 

STEP Cit3lJ>Join the two bignodes obtained above by an edge having 
a capacity equal to the maximum flow value computed in 
(i>. This edge is the first link of the cut-tree. The 
cut-tree also has two bignodes now. 

STEP C ivJ> Select any bignode in the cut-tree developed so far, 
which contains two or more vertices. Repeat steps 

by choosing a SOURCE and a SINK from the 
bignode selected. 

If no bignode can be selected in step t . all 

bignodes contain only one vertex then the process 
halts. The cut-tree developed at the end of the 
process is the required Gomory-Hu cut-tree. 

For computing the maximum flow between SOURCE and SINK 
use a simplified network derived in the following manner :- 
Suppose the bignode selected in step C ix>:> is deleted from the 
cut-tree developed so far. This results in a disconnection of 
the cut-tree into various components. All the vertices belonging 
to any of the bignodes in one component are condensed into a 
single node. The vertices belonging to the selected bignode are 
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not condensed and the edges between them are copied from the 
original graph. Edges are drawn between two condensed nodes or a 
condensed node and an original vertex. The capacity of an edge 
drawn between two condensed nodes is equal to the sum of the 
capacities of the edges in the original graph joining a vertex in 
one condensed node to a vertex in the other condensed node. 

Similarly, the capacity of an edge drawn between a condensed node 

and an original vertex is equal to the sum of the capacities of 
the edges in the original graph joining a vertex in the condensed 
node to the original vertex. (The rationale behind the process 
of condensing a subset of vertices into a single vertex for the 
purposes of simplifying the maximum flow computations is 

explained in Section 2.5) 

It IS seen that the above procedure for constructing a 
cut-tree creates a partition of the vertices in each iteration. 
By allowing a partition in step CttJ) only when the maximum flow 
value computed in step CiD is smaller than K, the cut-tree 

obtained at the end of the procedure will contain bignodes which 
are partitions such that .the maximum flow between any two 

vertices contained in a bignode is greater than K. This being 
the case, there is no need to go for a separate partitioning 

algorithm like the one described under Section 2.1. However, it 

should be kept in mind that such a short-cut procedure may 

require, for each bignode to be partitioned, maximum flow 
computations for more than one pair of vertices in the bignode. 
If there are n. vertices in the bignode. O(n^) maximum flow 

I 
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computations may be needed for resolving the problem of 
partitioning the i bignode. On the other hand, the complete 
cut-tree can be constructed by solving n-1 maximum flow problems. 
The short-cut procedure does, however, provide an alternative to 
the earlier scheme of partitioning, based on processing the 
complete cut-tree to get maximum flow values for all pairs of 
vertices, and it's implementation may be more efficient in some 

cases. The cut-tree generated at the end of the procedure with 

the above modification is called a partial cut-tree of edge 

connectivity K. The algorithm for partial cut-tree can also be 
used to generate the complete cut-tree by setting the parameter K 
= n-1. The algorithm for constructing a partial Gomory-Hu 
cut-tree of specified edge connectivity K, is derived by 
modifying the procedure given by Gomory and Hu C3Z] for 

developing a complete cut-tree. The algorithm for constructing 
the partial cut-tree is given below. 

ALGORITHM PART I AL_CUT_TREE 
DATA STRUCTURES 

C_TREE : Vertex adjacency list format for storing the 

cut- tree 

WGR : Vertex adjacency list format for storing the 

working graph. The work i ng . graph is derived 
from the input graph GC with some subsets of 
the vertices condensed to single vertices, 
wherever possible, for the purpose of 
Maximum-Flow computations 
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MATRIX_ : A n X n matrix to store Maximum-Flow values 

MAXFLOW between all pairs of vertices ( SDURCE-SINK 

pairs ) 

BIGNODE : A node of C_TREE which contains ail the 

vertices of the partition induced subgraph it 
denotes 

N_VERT : Number of vertices in WGR 

CONTINUE_ 

SPLITTING : A boolean variable to store a decision 

STEP 1 INITIALIZE C_TREE 

STEP 1 INITIALIZE WGR 

STEP 5 IF N_VERT > 1 THEN CONTI NUE_SPL ITT ING := true 
STEP 4 WHILE CONTINUE_SPLITTING DO 

STEP 4.1 SELECT SOURCE and SINK vertices 

IDENTIFY the BIGNODE in C_TREE which contains 
SOURCE and SINK vertices 

IF unable to maKe any selection -THEN 

CONTINUE_SPLITTING := false 
STEP 4.2 IF CONTINUE_SPLITTING THEN 

STEP 4.2.1 CONDENSE all disjoint subsets of 

vertices not in BIGNOOE into single 
vertices in WGR, wherever possible 
STEP 4.2.2 DO MAXIMUM-FLOW computation on WGR 

between SOURCE and SINK to get MAXFLOW 
value 

STEP 4.2.3 UPDATE values in MATRIX_MAXFLOW for the 

SOURCE-SINK pair 



STEP 4.2.4 


IF MAXFLOW < K THEM MODIFY C TREE i.e. 


create a partition induced subgraph in 
the form of additional bignodes in 
C_TREE 

STEP 4.2.5 RE-INITIALIZE WGR 

3.2.3 Condensing a subset, of vert-ices into a single vertex 

Maximum flow computations required to be done in the 
ALGORITHM PART I AL_CUT_TREE , typicafiy consume a time bounded by a 
polynomial function of the number of vertices n and the number of 
edges m in the graph. For example, Karzanov's maximum flow 
algorithm runs in 0<n^) time. The maximum flow computations can 
be done on a simplified graph with lesser number of vertices, 
called the working graph.. This will result in a reduction in the 
number of computations to be made. The working graph WGR is 
derived from the original graph GC by utilising a process of 
condensing a subset of vertices into a single vertex. An example 
is illustrated in Figure 5.6. 

By condensing 4, 5 and 6 into a single vertex the edge 
capacities of <4,5>, (4,6) and (5,6) have been raised from 1 to 
infinity. Since these edges are not in the 2-5 minimum cut, any 
increase in their capacities does not affect the maximum flow 
value between the vertices 2 and 5. Therefore, so far as 
calculating the maximum flow value between the vertices 2 and 5 
IS concerned, the vertices 4, 5 and 6 may be condensed into a 





Condense vertices 4.5 and 6 


(b) Simplified graph WGR 


NOTE : - The maximum flow between 2 and 3 is 3 units and it is the same 
in both GC and WGR. 


FIGURE 3.6 Condensing a subset of vertices into a single vertex. 



single vertex and the computations may be done on the graph WGR. 


There are a couple of conditions under which the 
process of condensation is allowed and they are stated below in 
the form of Lemma <5.1> and Lemma <5.2> given in CHu,*]. 

[Lemma 3.11 Let (X, X) be a minimum cut separating v^ e X 
and some other vertex, and let v and v. be any two vertices 
contained in X. Then there exists a minimum cut (Z, Z) 
separating v and v such that <Z. Z) and (X, X> do not cross 
each other. ■ 

(Refer to Figure 5.7 for details of the cuts.) 

This means that if (X, X) is a minimum cut, the subset 
X can be condensed into a single vertex for computing the maximum 
flow between any two vertices in X. 

[Lemma 3.21 Let (X. X> be a minimum cut separating v^ and 
some other vertex, and let v, be any vertex that belongs to X. 
Then there exists a minimum cut (Z, Z) separating v^ and v^ such 
that (X, X) and (Z, Z) do not cross each other. ■ 

(Refer to Figure 5.7 for details of the cuts.) 

This means that if (X. X) is a minimum cut separating 
V. and some other vertex and the value of the maximum flow 
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between i and k is to be determined, where v, is any vertex in X, 
then the subset X can be condensed into a single vertex. The 
STEP 4.2.1 in ALGORITHM PARTI AL_CUT_TREE for condensation is 
carried out based on the implications of Lemma (5.1) and Lemma 
(3.2) to get a simplified working graph WGR. The maximum flow 
computations are carried out on WGR to effect some savings in the 
amount of computation required. 


3.2.4 Algorithm for Maxiimim Flow computation 

Any of the standard routines, surveyed below, can used 
for computing maximum flow between a source vertex and a sink 
vertex. The basic steps of a maximum flow algorithm are 

(i) Based on the current flow, the residual network is 
constructed. Each arc in the residual network has a 
non~zero capacity, called residual capacity = capacity 
of the arc in the original network - flow along the 
arc . 

(ii) A flow augmenting path (path along which flow can be 
sent) between the source vertex and the sink vertex in 
the residual network is selected and flow is sent along 
this path. 


For arbitrary selection of the flow augmenting path. 
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the running time of the above algorithm can go upto 0<nmU) where 

n is the number of vertices, m is the number of edges in the 

graph and U is the upper bound on the arc-flow capacities. The 

first major improvement in algorithmic efficiency came from 

Edmonds and Karp who selected the shortest flow augmenting path 

(i.e. the path with the least number of arcs) every time in step 

(ii). The running time of this algorithm is bounded by OCnm >. 

There are many algorithms which improve upon this bound. Most of 

these are based on the algorithm given by Dinic which uses the 

concept of layered networks. Dinic's algorithm runs in 0(n m) 

time. A major improvement upon this was given by Malhotra &i a.1 . 

Their MPM algorithm has a running time bounded by 0<n^). A 

fundamentally different approach using the method of pre-flows is 

used by Karzanov’s algorithm which is also 0(n^> in computational 

complexity. The technique of capacity scaling adopted by Tarjan 

et al has further pushed the bound which is practically 

2 

achievable to OCnm log<n /m)). 

Keeping in view that most of the time arc-flow 
capacities in the graph GC are small numbers, the basic flow 
augmenting path algorithm can be expected to do well from a 
computational point of view. The implementation of this 
algorithm is also considerably easier. 
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3*3 Merging Algoriihm for Colored Subgraphs 

Any generalized algorithm for merging two colored 
subgraphs can be used to merge more than two colored subgraphs by 
an iterative process. The generalized algorithm for merging two 
colored subgraphs GR_A and GR_B takes them as input along with a 
list L of edges which join a vertex in GR_A to a vertex in GR_B 

in the original graph GC. The output is a colored subgraph 

GR_A_B induced by the vertices of GR_A and GR_B . In general, 
GR_A_B uses more colors than GR_A and GR_B . While the ordering 
of the edges in the input list L may be arbitrary, the number of 
colors added during the merging process depends upon the process 
itself and the ordering of the vertices in the list L. 

The algorithm is structured such that each edge in the 
list L has a head which is a vertex of GR_A and a tail which is a 

vertex of GR_B . Each edge in the list L ,is examined once in the 

order in which it appears. Initially, the coloring for each 
vertex in GR_A and GR_B is unconstrained i .e . it is free to be 
changed during the process of merging. One of the following 
steps IS taken upon examining each edge in the list L 

(i) If COLORChead) COLOR! ta i 1 ) then the coloring for the 
head and the tail are constrained to their respective 
colors. 

(ii) If COLOR<head) = COLOR(tail) and coloring of both head 
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and tail are constrained to their respective colors 
then an additional color is used for the head ( or 
alternatively the tail }. 

(iii>If COLOR<head) = COLOi?(tail) and coloring of either 
head or tail is constrained, say that of head, but not 
both or coloring of neither is constrained then the 
unconstrained colors in the subgraph GR_B ( which 
contains the tail ) can be re-arranged to get a 
re-coloring of the tail which does not conflict with 
the coloring of the head. 

ALGORITHM MERGING 

STEP 1 : (Initialization) 

FOR I := 1 TO N DO CONS_ACI] := FALSE 

CONS_BCI] := FALSE 

/* Color of vertex I in GR_A and GR_B is 

unconstrained initially */ 

TEMPEDGE := PO I NTER_TO_TOP_OF_L I ST ( L) 

/* L is the list of edges which join a a vertex in 
GR_A to a vertex in GR_B in the original graph 
GC , HEAD being in GR_A and TAIL in GR_B */ 
CQLOURS_A := NUMBER_OF_COLQRS_IN_GR_A 
STEP 1 : (Body of the algorithm) 

IF TEMPEDGE HIL /* List L is not empty */ 

/* One of the following four cases occurs */ 

CASE 1 ( COLOR < HEAD ( TEMPEDGE )_ I N_GR_A) 
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COLOR < HEAD <TEMPEDGE)__IN_GR_A) > 
CONS_ACHEAD(TEMPEDGE) ] := TRUE 
CONS_BCTAIL(TEMPEDGE) ] := TRUE 
/* Both head and tail are constrained to their 
respective colors */ 

CASE 2 ( COLOR ( HEAD (TEMPEDGE>_IN_GR_A> = 

COLOR<HEAD(TEMPEDGE)_IN_GR_A) > AMD 
CONS_ACHEAD<TEMPEDGE) ] = TRUE AMD 
C0NS_BCTAIL(TEMPEDGE) 3 = TRUE 
COLOR(HEAD<TEMPEDGE)_IN_GR_A> := C0LDURS_A+1 
/* An additional color is used */ 

CASE 3 < COLOR(HEAD(TEMPEDGE)_IN_GR_A> = 

COLOR ( HEAD (TEMPEDGE)_IN_GR_A> ) AMD 
CONS_A[HEAD(TEMPEDGE) 3 := TRUE AMD 
CONS_B[TAIL(TEMPEDGE) 3 = FALSE 
CALL RENAME_COLOR(CONS_B.GR_B> 

CASE 4 ( COLOR ( HEAD (TEMPEDGE)_IN_GR_A) = 

COLOR(HEAD(TEMPEDGE)_IN_GR_A> ) AMD 
C0NS_ACTAIL<TEMPEDGE)3 = FALSE 
CALL RENAME_COLOR(CONS_A.GR_A) 

TEMPEDGE := PO INTER_TO_NEXT_OF_TEMPEDGE 
GOTO / 2 

/* Next edge in the list L is examined */ 


PROCEDURE RENAME_COLOR(CONS,GR) 
STEP 1 : (Initialization) 


FOR I := 1 TO N DO CCI3 := 0 
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FOR I := 1 ro N DO CCCOLOR( VERETX_I_IN_GR) ] := 
COLOR (VERETX_I_IN_GR) 

FOR I ;= 1 TO N DO IF C0NSCI3 THEN 
CCCOLOR<VERETX_I_IN_GR) ] := 0 

STEP Z : (Rearrange colors if possible else use additional 

color) 

IF CCI ] 0 for more than one value of I THEN 

Rearrange in a cyclical order the values of 
CCI] amongst those I for which CCI] ^ 0 
COLOR (VERTEX_IN_GR> := CCI] 

/* Rearrange colors of vertices of GR */ 

ELSE Use an additional color for HEAD in GR_A as 
in CASE 2 of the main body of the algorithm. 

3. 4 Evaluation of Inter modulation Interference 

As explained in Chapter II. evaluation of the order of 
intermodulation interference is to be carried out at each site in 
the communication network for a given frequency plan. The 
evaluation problem was formulated as a nonlinear integer program 
(2.9) 

n 

M i n Q = E 
i = 

n 

Subject to 2 

il = 

n 

E 

i = 


F. X. < GB; 

1 i 

F. X. > -GB; 

i i 
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3 at least one i such that 


X . = 1 ; 

t 


are integers and free; 


where 1 < i < n. 


To declare the set of frequencies F,. F, F , which are 

1 c n 

active at a site in the network, to be intermodulat ion-free upto 

opt 

order Q < or less >, the optimal value of Q i.&. Q obtained by 

o o 

solving (.) should be greater than Q+1 . As there is no real 

• opt opt 

purpose in determining Q other than finding out whether Q > 

o o 

Q+1 or not, a bound on the value of may be placed through an 


additional constraint 


Q < 
o 


(3.12) 


If there is no feasible solution to <2. 9) along with the above 

opt opt 

constraint then it means that Q > Q+1, else Q < Q+1. 

o o 


The above nonlinear integer program can be linearised 
easily by introducing fresh variables X^ and X^ such that 

X. = X'^ - xT (3.15) 


'he resulting linear integer program can be written as 


Min E (X. + X ) 

1 = 1 


(5.14) 


Subject to 


r F,(X.-X.)< GB; 

, 1 . ^ 

1 = 1 


(3.13) 


r F . < X - X . ) > -GB; 

, i I i 

1 = 1 


(3.16) 


3 at least one i such that X^ ^ ~ ^ ' 


(3.17) 
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, X^ > 0 and integers, where 1 < i ^ n. 


The optimal solution of the LP-re laxat i on of the above problem is 

one of the X^ , xT = 1 ; 

+ ^ • 

X = -=1- ; and <5.18) 

J 



rest of the X^ , X, = 0, 
is the highest frequency among the n-1 
and 


where k i , j and, 
frequencies excluding 


LP-opt F. 

Q = 1 + <3. 19) 

o F . 

J 


LP-opt 

In most practical cases of intermodulation problems, is 

opt 

substantially smaller than the optimal integer solution 
Therefore any attempt at a solution by using a branch and bound 
procedure runs into the difficulty of obtaining a good lower 
bound based on the LP— re laxat i on . So backtrack algorithms are 
used as a solution strategy. 


In the primal backtrack algorithm, an initial solution 
is found and then improved solutions are sought. The procedure 
terminates when no solution with a lower Q^-value can be found. 

The dual backtrack algorithm which has been used to 

solve the problem (2.9) along with the constraint <3.12), starts 

with Q = 1 and tries to find a feasible solution to the problem 
o 

for the fixed value of Q^. failing which it increments by 1 
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and repeats the solution procedure until it fjnds a solution or 

= Q + 1 . If it has found no solution upto the stage Q = Q+1 

opt 

then it exits with the inference that Q > Q+1. The first 

o 

solution, if found, is also the optimal since it is the dual 
problem which is being solved. The algorithm runs as follows 

ALGORITHM DUAL_BACKTRACK 

STEP 1 : < Initialization > 

GOAL : = 1 ; 

Xj := - GOAL; 

X,. X^ := 0; 

L y n 

LEVEL := 2; 

ONE := 1; C ONE indicates the number of s 

fixed at the value ±1 } 

GOTO / 1 . 

STEP 2 : ( Optimality Check > 

n 

IF ( ABS < E F. X. ) < GB ) AMD 

i =1 

( ONE > 1 > THEM STOP 
ELSE GOTO f 2. 

STEP 5 : ( Construct next solution candidate > 

LEVEL 

IF 1^1 “ GOAL THEN GOTO / 2c. 

i = 1 ^ 

STEP 2a : ( Forward step ) 

IF LEVEL = n then := - 

GOTO / 1 

ELSE := + 1 



MODIFY ONE 


STEP 


STEP 


2b 


LEVEL := LEVEL + 1 

LEVEL - 1 

^LEVEL ^ 

1 = 1 

MODIFY ONE 


GOAL 


GOTO / 2b. 

( Implicit enumeration tests > 


LEVEL - 1 

IF ( E f.X. 

. , VI. 


'^LEVEL \eVEL 


2c 


< i > 


OR 

LEVEL - 1 

^ + ^EVEL \eVEL > 

I. = 1 

THEM GOTO /2c 
ELSE GOTO / 1 . 

( Backward step > 

IF LEVEL = 1 THEM GOTO J 5 

1 

+ 1 . 


ELSE : = 0 

MODIFY ONE 
LEVEL := LEVEL 


\evel \evel 


MODIFY ONE 
LEVEL 

( 1 1 > IE E I I = 

i = 1 

THEM GOTO t 2c ( i ) 

ELSE LEVEL := LEVEL +• 1 
LEVEL - 1 

^LEVEL 


1 


X . - GOAL 



MODIFY ONE 


GOTO / 2b. 

STEP 5 : (No optimal solution with Q 

o 

GOAL := GOAL + 1 

IF GOAL < Q THEM GOTO / 1 . 


= GOAL } 



CHAPTER IV 


IMPLEMENTATION AND RESULTS 

The solution niethodologies , which were outlined in 
chapter III, for the various sub-problems of the Frequency 
Planning Problem (FPP) are used to develop computer programs 
which will facilitate quick solutions. In this chapter, the 
organization of the various computer programs is explained along 
with the functional aspects of each program. The data handling 
done by the various programs, using disk file structures, is also 
explained for the purposes , of understanding the input 
requirements and interpreting the output results. This is 
followed by an example of frequency planning using the system of 
computer programs and evaluation and interpretation of the 
results obtained. Directions for further work and improvements 
are explored as a part of the conclusion. 

4# 1 Over-view of ihe Frequency Planning System 

The process of frequency plannihg, as envisaged in this 
thesis, IS depicted in the flow diagram in Figure 4.1. . Three 
data files are assumed to exist at the beginning of the process - 
SITE.DAT, GR.DAT and COL. DAT. Information regarding incidence of 
the links at various sites is maintained in SITE.DAT. The graph 




FIGURE 4.1 Flow diagram showing the frequency planning process 






76 


to be colored- GC - is stored in an adjacency list format in 
GR.DAT. The data file COL. OAT stores the one-to-one mapping of 
colors onto the set of frequencies available for use by ,a 
communication link. The planning process commences by reading 
into the memory SITE_GR from SITE.DAT and the graph GC from 
GR.DAT, both the data structures being arrays of singly linked 
lists as defined in PASCAL. Next, the FREQUENCY PLANNER module 
utilizes a graph coloring strategy to color the graph GC. The 
coloring strategy could be to partition the graph GC. color all 
the components and merge the colorings of the components to get a 
coloring of GC. Alternatively, the entire graph GC could be 
colored by a heuristic coloring algorithm. The frequency 
assignment for the various links is stored in FREQ.DAT. In the 
next stage, the user- i nteract i ve FREQUENCY PLAN EVALUATION module 
is used iteratively to finalize a frequency plan. In this 
module, the given frequency assignment is evaluated for 
occurrence of a significant order of intermodulation at any site 
in SITE_GR and based on the results, subsequent user-modification 
of frequency assignment to the links is permitted. A frequency 
plan modified by the user must undergo the evaluation described 
above. When no modification remains to be done, the frequency 
plan is stored in PLAN.DAT. 

4. 2 Organization of Conq>utor Programs and Files 

The computer programs, for implementing the solution 
methodologies for the various sub-problems of FPP, have been 
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coded in PASCAL using the software package Turbo-Pascal 5.0 and 
can be run on an IBM-compatible PC. The facility for defining 
units has also been exploited to enhance the modular ana 
structured approach to programming, in general. A similar 
facility exists in HP-Pascal for defining modules and has also 
been utilized to enable the programs to run on HP super-mini 
computer systems. 

There are two unit files which contain the various 
connmonly used procedures. These are the unit LINK. PAS and the 
unit FP.PAS. The unit LINK. PAS contains the declarations 

regarding the data structures common across all the procedures, 
notably that of the GRAPH data. stTXicfure. It also contains a 
procedvjre R_GRAPH which reads the formatted input from a text 
file and stores it in a GRAPH . This is nothing but the input 
undirected network in adjacency list format. 

The unit FP.PAS has comparatively more number of 
procedures. These are :- 

( 1 ) procedure MAX_FLOW : this solves the maximum flow 

problem on the specified network WGR, between the nodes 
SOURCE and SINK and gets, the' maximum flow value as well 
as the minimum cut partition of the nodes of V/GR. 
There are four sub-procedures contained in HAX_FLOW . 
BFS DISTJLABEL performs a breadth-first-search of the 
working network WGR starting from node SINK to obtain 
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dista.nc8 labels in the array DIST. The sv.b>—pi'ocec£uT& 
ADM1SSIBLE_ARC finds out an arc, if it exists, from the 
piesent node which is admissible in the augmenting path 
from SOURCE to SINK. RETREAT backtracks from the 
current node if none of the arcs incident on it is 
admissible. AUGMENT increases the flow along the 
augmenting path identified between SOURCE and SINK. 


( 1 i ) procedxire CUT_TREE uses the procedxxre MAX_FLOW to 
develop a partial Gomory-Hu cut-tree of specified edge 
connectivity K. The cut-tree is stored in a GRAPH 
called C_TREE , in adjacency list format. There are 
many sub-procedures in CUTJTREE. INIT_C_TREE 

initializes the adjacency list format of C_TREE into a 
single bignode. The SOURCE-SINK node pair is chosen by 
SELECT _SOURCE_S INK from the same component (bignode) of 
C_TREE. INIT_WGR initializes the network WGR into a 
null network. CONDENSE carries out the process of 
condensing all nodes in one component of the input 
graph GC as indicated by the cut-tree CJTREE and 
constructs a simplified network WGR. MODIFY CUT TREE 


modifies the CJTREE after the solution of the nmximum 
flow problem between SOURCE and SINK yields an 


acceptable partition of the nodes. 


WRITE GRAPH 


writes the CJTREE developed at any stage to the output 


file in adjacency list format. 
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< i i i IprocedT-ire GR_COL ; colors the input graph using the 
sequential coloring algorithm for SL as well as LF 
ordering of the vertices. It has two sub-procedures 
called ORDERING and SEQCOLORING. ORDERING generates a 
SL/LF ordering of the vertices while SEQCOLORING colors 
the vertices sequentially in the order given by 
ORDERING . 


<iv> pTocedure WR_GRAPH : writes the colored graph to the 
pre-defined output file. 

<v) procedure COL_FREQ : reads the frequencies 

corresponding to the colors used by the procedure 
GR_COL , from the file COLOR.DAT. 

(vi > procedure WR_FREQ_PLAN : writes the current frequency 
plan to the file FREQ.DAT. 

The program COLOR. PAS uses the procedures R_GRAPH , 
GR_COL , COL_FREQ and WR_FREQ_PLAN to generate the first frequency 
plan in the file FREQ.DAT. As there are two possible orderings. 
SL and LF , of the vertices in GR_COL ^ two first frequency plans 
can be generated. The program EVAL.PAS evaluates the current 
frequency plan for critical cases of intermodulation interference 
and allows the user to modify it in an interactive manner to 
reduce the order of intermodulation. It uses the following 


procedures 



W 


(i) procedure READ_SITE_GR ; reads the corrHinun i cat i on 

network as a link-site incidence list from SITE.DAT 

into SITEJ3R. 

(li) procedux-e READ_LIHK_FR ■. reads link frequencies from 
FREQ.DAT which is the first frequency plan, into array 
LIMK_FREQ. 

{ i i \ ) procedux-e EVAL_IMTERMOD : evaluates the frequency 

assignment by checking for intermodulation of order 

less than MAXORDER at each site in SITE_CR. A dual 
backtrack algorithm is used for this purpose. 


(iv> procedux-e EV AL_SEPARATIOM : releases a list of 
frequencies in an array AV _FR£1Q , which can be used for 
assignment to a specified link without violating 
the co-channel or adjacent-channel constraints. This 
list IS used to prompt the user in modifying the 
frequency plan. 


(v) procedure FPJiODIFY : allows the user to change the 
frequency assignment to the specified link. 

The program PARTCOL.PAS uses the procedure CUT_TREE to 
generate a partition of the input graph GC . Next, it invokes the 
sub-procedures in GR_COL to color each partition individually. 
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Then pTocedure RE!V_COL is applied which, given two colored 

^ ""9 ^'3-Ph s CR. A and CR R, and a list L of edges which join a 

vertex in GR_A to a vertex in GR_B , gives a revised coloring of 
the vertices of the sub-graph formed by the merger of GR_>4 and 

Repeated applications of RE!V COL. are used to completely 

color GC and thus produce an al t&i'ncctix>& first freqxi&ricy plan. 
This frequency plan is written into FREQ.DAT and subjected to the 
frequency assignment evaluation program EVAL . PAS . 

The listing of files is as given below 


PROGRAM FILES 

DATA FILES 

LINK. PAS 

COLOR.DAT 

FP.PAS 

SITE.DAT 

COLOR. PAS 

GR . DAT 

PARTCOL.PAS 

FREQ.DAT 

EVAL.PAS 

PLAN.DAT 


REPO.DAT 


4. 3 A Solved Example using the Frequency Planning System 

The system of computer programs described earlier has 
been used to develop a frequency plan for a portion of the Indian 
Railways MW network. Table 4.1 shows the contents of the file 
COLOR.DAT. These are the frequency resources onto which the 
colors map uniquely. The portion of the IR MW network taken up 
for frequency planning is shown in Figure 4.2. The data for the 



TABLE 4..1 Frequency Resources available for use 


Co 1 or 

No . 

Corresponding Microwave 

Lower Band 

Frequency' (MHz) 

Upper Band 

1 

7128 

7289 

2 

7155 

7296 

5 

7142 

7505 

4 

7149 

7510 

5 

7156 

7517 

6 

7165 

7524 

7 

7170 

7551 

8 

7177 

7558 

9 

7184 

7545 

10 

7191 

7552 

1 1 

7198 

7559 

12 

7205 

7 566 

15 

7212 

7575 

14 

7219 

7580 

15 

7226 

7587 

16 

7255 

7594 

17 

7240 

7401 

18 

7247 

7408 

19 

7254 

7415 

20 

7261 

7422 







FIGURE 4.2 The example network stored in SITE_GR 



FIGURE 4.3 The graph GC to be colored 





i fie 


iinK-site incidence representation of the 
network of Figure 4.1 


Site No . 

Link No.s incident at this site 

1 

1 2 5 4 5 

2 

16 7 8 

5 

6 9 

4 

9 10 11 12 15 


10 14 15 

6 

14 16 17 18 

7 

2 7 19 20 

8 

19 21 

9 

21 22 

10 

20 25 

11 

25 24 

12 

5 24 

15 

4 25 

14 

25 

15 

8 26 

16 

26 27 

17 

27 

18 

16 28 

19 

28 

20 

1129 

22 

50 52 

25 

12 15 17 51 55 

24 

18 52 55 

25 

15 22 

26 

5 



TABLE 4-3 Adjacency list representation of graph GC to be 


colored 


Link No. 

Adjacent Links 

i .e 

. interfering links 

1 

2 

4 

6 

7 

8 


2 

1 

3 

6 

7 

8 

19 20 

5 

2 

5 

20 

24 



4 

1 

6 

8 

25 

26 


5 

3 






6 

1 

2 

4 

7 

8 

9 15 26 

7 

1 

2 

6 

8 

19 

20 21 

8 

1 

2 

4 

6 

7 

26 

9 

6 

10 

11 

15 



10 

9 

11 

12 

14 

15 

29 

1 1 

9 

10 

12 

15 

29 


12 

10 

11 

14 

15 

29 

51 

15 

6 

9 

1 1 

22 



14 

10 

12 

1 5 

16 

17 

18 

15 

10 

12 

14 

17 

18 

29 51 

16 

14 

28 





17 

, 14 

15 

18 

51 

55 


18 

14 

15 

17 

51 

55 


19 

2 

7 

20 

21 



20 

2 

5 

7 

19 

25 

24 
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TABLE 4.3 contd. 


21 

7 19 22 

22 

15 21 

25 

20 24 

24 

5 20 25 

25 

4 26 27 

26 

4 6 8 25 27 

27 

25 26 

28 

16 

29 

10 11 12 15 50 51 

50 

29 51 52 55 

51 

12 15 17 18 29 50 52 55 

52 

50 51 55 

55 

17 18 50 51 52 



TABLE 4-A The First Frequency Plan (LF ordering) 



7191 

7163 

7149 


4 

3 

6 

7 

8 
9 

10 

11 


15 

16 

17 

18 
19 


7165 

71Z8 

7155 

7149 

7177 

7165 

7184 

7155 
7170 
7149 

7156 
7142 
7155 
7184 
7170 
7177 


20 


7135 




TABLE 4.5 A Report on Intermodulation Evaluation of Frequency 
Plan given in Table 4.4 


SITE 1 Order of intermodulation = 2 

FREQUENCIES IN INTERMODULATION : 7165 7165 

MULTIPLIER XCJ] : -1 1 

LINK NO. : 4 2 

No. of Iterations for solution = 144 


SITE 


2 

Order of 

i ntermodu lat ion is 

high I 

SITE 

// 

5 

Order of 

intermodulation is 

high ! 

SITE 

// 

4 

Order 

of 1 ntermodu lat ion 

= 5 


FREQUENCIES IN INTERMODULATION : 7149 7170 7184 

716 5 

MULTIPLIER XCJ3 : 1-11 

-1 

LINK NO. : 15 12 10 

9 


No. of Iterations for solution = 651 


SITE 

// 

5 

Orde r 

of 

1 n t e rmodu 1 at i on 

is high ! 

SITE 

// 

6 

Order 

of 

i ntermodu lat i on 

is high ! 


SITE 7 Order of i ntermodu lat i on 

FREQUENCIES IN INTERMODULATION : 
MULTIPLIER XCJ] = 

LINK NO. : 

No. of Iterations for solution = 227 


= 5 

7155 7149 7165 

-1 2 -1 

20 7 2 


SITE // 8 


Order of intermodulation is high ! 


SITE 

// 

9 

Order 

of 

i nt ermodu lat ion 

i s 

h i gh 

J 

SITE 

// 

10 

Order 

of 

intermoduiat ion 

i s 

h i gh 

J 

SITE 

# 

11 

Order 

of 

i nt ermodu lat i on 

i s 

h i gh 

1 

SITE 

4? 

12 

Order 

of 

i nt ermodu lat i on 

i s 

h i gh 

t 

SITE 

// 

15 

Order 

of 

i nt ermodu lat i on 

i s 

h i gh 

; 

SITE 

// 

14 

Order 

of 

i nt ermodu iat ion 

i s 

h i gh 

J 

SITE 

// 

15 

Order 

of 

i ntermodu lat i on 

i s 

h i gh 

t 

SITE 

// 

16 

Order 

of 

1 n termodu lat i on 

i s 

h i gh 

j 

SITE 

77 

17 

Order 

of 

i ntermodu iat i on 

i s 

h i gh 

I 

SITE 

// 

18 

Order 

of 

i ntermodu lat i on 

l s 

h i gh 

J 

SITE 

// 

19 

Order 

of 

i ntermodu iat i on 

1 s 

h 1 gh 

; 

SITE 


20 

Order 

of 

i ntermodu lat i on 

i s 

h i gh 

1 

SITE 

// 

21 

Orde r 

of 

intermoduiat i on 

l s 

h i gh 

\ 

SITE 

// 

22 

Orde r 

of 

intermoduiat i on 

1 s 

high 

1 

SITE 

7^ 

25 

Order 

of 1 ntermodu lat ion 

= 2 



FREQUENCIES IN INTERMODULATION : 7142 7142 

MULTIPLIER XCJ] : "1 i 

LINK NO. : 55 15 

No. of Iterations for solution = 67 



SITE // 24 


Order of intermoduiat ion 


FREQUENCIES IN INTERMODULATIQN : 
MULTIPLIER XCJ3 : 

LINK NO. ; 

No. of Iterations for solution = 145 


SITE // 25 


SITE # 26 


Order of intermodulation is high 


Order of intermodulation is high ! 


TABLE 4-.6 The Final Frequency Plan <LF ordering) 

(after 25 iterations) 


L 1 nk 

Frequency 

No . 

(MHz) 

1 

7165 

2 

7205 

5 

7177 

4 

7212 

5 

7156 

6 

7128 

7 

7149 

8 

7177 

9 

7165 

10 

7184 






TABLE 4-.6 contd. 


11 

7155 

12 

7170 

1 ? 

7177 

14 

7156 

15 

7155 

16 

7155 

17 

7184 

18 

7170 

19 

7177 

20 

7155 

21 

7155 

22 

7165 

25 

7149 

24 

7165 

25 

7155 

26 

7149 

27 

7165 

28 

7149 

29 

7156 

50 

7170 

51 

7226 

52 

7156 

55 

7198 




9 ^ 


equivalent SITE_GR (stored in SITE. DAT) is given in Table 4.2. 
Based on interference ca 1 cu lat i ons . the graph GC to be colored is 
drawn up in Figure 4.5 and it‘s adjacency list representation 
(stored in GR.DAT) is given in Table 4.5. 



The first frequency plans 

produced 

(for LF ordering of 

1 i nks ) i s 

displayed in 

Table 4 

.4. 

A 

report on the 

1 ntermodulat ion evaluation 

appears 

i n 

Tab 1 e 

4.5. The final 

frequency 

plan, produced 

after 

23 

number 

of iterations of 


user-modifications, is presented in Table 4.6. 

4.4 Conclusion 

The present thesis modelled the problem of frequency 

planning as a graph coloring problem. To color the input graph, 

♦ 

the standard sequential algorithm is used which is an approximate 
one. The comparative performances of graph coloring algorithms 
have been evaluated in literature ( [26], [35]). Consequently, a 

better algorithm like the RLF [55] may be implemented to generate 
a better first frequency plan. In fact, one of the advantages of 
modelling the FPP as a graph coloring problem is that any 
coloring algorithm developed is directly usable for solving the 
FPP. 


The alternative strategy of partitioning the input 
graph, coloring the partitions and merging them, involves the 
development of a partial cut-tree. The maximum flow problem 



solved therein used the Shortest Augmenting Path Algorithm. 

\ 

Since the arc capacities are generally small numbers in this 
case, the time spent on finding the shortest augmenting path may 
be saved by using any augmenting path instead. The extra time 
spent on the increased number of flow augmentations may be offset 
by the savings. In addition, other partitioning techniques like 
those based on removal of vertex sets or cliques, may be taken up 
for consideration. 

The i nt ermodu lat i on evaluation is carried out 
iteratively by the user to modify the first frequency plan and 
eliminate intermodulation interference. This is because all the 
possible frequency assignments have to be implicitly or 
explicitly subject to i nt ermodu lat i on evaluation. The evaluation 
problem is itself a nonlinear integer program and it’s solution 
incurs considerable computational cost. Hence, judicious 
intervention on the part of the user helps in modifying the first 
frequency plan in a small number of iterations. A method of 
ordering the links and frequencies may be sought before 
generating the first frequency plan such that the number of cases 
of intermodulation violation is minimal. It may be noted that 
the sequential coloring algorithms do order the links in a 
particular fashion -(SL/LF) before coloring them. 

A branch and bound procedure for generating an 
acceptable frequency plan is feasible for small networks and it 
may be tried out. Such a scheme works by taking a forward step or 
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a backward step at 

each stage. At any 

stage 

there exists 

a 

part lai 

ass i gnment 

of frequencies to some 

of the 

1 i nks . 

In 

a 

f orward 

step, one 

more link is assigned a 

frequency; while 

i n 

a 


backward step, the previous assignment is revoked. After each 
forward step, the partial assignment is checked for the 
interference level in each link (Formulae, like those given in 
Appendix I, will be useful). An upper bound on the interference 
level is provided by assuming that ail the unassigned links have 
the same frequency (maximum possible interference). After each 
backward step, all the infeasible partial assignments are pruned 
as they can never lead to a feasible frequency plan. The 
ordering of links and frequencies can be very crucial to the 
running time of such a procedure. 
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APPENDIX I 


1 A technical overview of the operation of a link in 
a commxinica'tion network 

The scheme of operation of a po i nt -to-po i nt link in a 
communication network is illustrated for the case of a Miarawa\>& 
system, (See Figure I.l) The process of a signal being relayed 
from Site 1 to Site 2 involves the following steps 

(i) A si^al to be relayed is generated by the concerned 
equ 1 pmen t . 

<ii) The carrier xixive (it has a frequency in the microwave 
spectrum range, which is very large compared to the signal 
wave frequency) is generated using a special device called 
oscillator which could be a vacuum tube device like the 
klystron , 

(iii>The mjodxilator combines the signal wave with the carrier wave 
and feeds the output to the transmitter after filtering out 
spurious emissions at frequencies outside the neighbourhood 
of the modulated wave frequency. The effect of filtering is 
to reduce the effect of adjacent channel interference 
(described under. Sec t i on 5 ) but the filtering may not be 


adequate always. 



downlink 


uplink 



FIGURE 1.1 Microwave communication link 
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(iv) The transmitter boosts the power of it’s input wave (signal 

+ carrier) upto the required levels and feeds it to the 

antenna via the \j>aveguide . 

<v> The antenna in turn, by the virtue of it's geometry, beams 
the input wave in the direction of the receiving antenna. 

(vi> The receiving antenna feeds the received wave to the 

demodulator via the waveguide. 

(vii)The demodulator separates the signal wave from the carrier 
wave using a process similar but inverse to that of 

modulation. It also employs filtering to remove interfering 
signals at unwanted frequencies. 

Wherever, the text of the thesis contains a mention of 
frequency planning, the frequencies being referred to are 
actually those of the carrier waves which fall in the microwave 
range of the spectrum. Depending upon the type of service 

provided by the communication network, it can avail of the 
frequency bands (and channels thereof) allocated to that type of 
telecommunication services at the international and national 
levels. An additional consideration would be the previous 
frequency assignments made in the same frequency band, if any, 
for an existing network providing similar telecommunication 
services. So the frequency planner for the communication network 
ends up with a small set of discrete frequency channels out of 
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the entire usable spectrum. This set of frequencies constitutes 
the electromagnetic resources shared by all the links of the 
communication network. 


A sirnultaneoTjLS t-wo-VKzy communication between the two 
sites on a link, called a duplex link, will require two different 
frequencies for the up-iink and the down-link respectively, 
because the signal paths for up-link and down-link are the same 
and the two signals can potentially interfere with each other. A 
radio telephone link is an example of a duplex link. A 
half— duplex link, on the other hand, allows two-way communication 
on a single link frequency but simultaneous operation of the 
up-link and the down-link is not permissible, as for example in 
the case of wireless sets. The simplex link is a one-way 
communication only, using a single link frequency, as in radio 
broadcasting networks. 


If a site is linked to more than one site, then a 
common antenna and waveguide might be used in the transmitter / 
receiver system to meet the requirements of some or all the links 
incident at this site. The antenna geometry is usually paraboltc 
with the signal wave being fed at it's focus, either directly 
from a horn or via a reflector . The signal wave at the focus of 
the parabolic antenna gets transmitted as a parallel beam after 
reflection from the parabolic surface of the antenna. The 

directionality lent by the antenna to the wave beams restricts 
their spherical dispersion in the atmospheric media and reduces 
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the spillover of the ( i nterf ering) signals into *wrong’ 
receivers. This causes an attenuation of the signal in the 
directions different from the direction in which the signal wave 
is beamed and is called angular attenuation. In such a case, a 
greater re-use of frequencies may be possible and this leads to 
increased utilization of the frequency spectrum. 

Polox-ixatton of a wave restricts the wave motion to a 
plane. (See Figure 1.2) If two non-po lar i zed waves, incident at 


a receiver, are 

potent ial ly 

interfering 

with 

each 

other 

then 

cross-polarization 

of the 

two waves 

(the 

two 

waves 

are 

plane- polarized 

i n 

or thogona 1 

planes) < 

wi 1 I 

reduce 

the 

mutua 1 


interference potential. So frequency re-use can be increased by 
resorting to this technique. 


2: Transmission Losses during the propagation of a radio wave 


In any communication system, the transmitted wave 
suffers a loss in power as it propagates away from the 
transmitter and reaches the receiver. To quantify the amount of 
loss taking place, a few definitions are introduced below 


A linear t ime- i nvar lant communication system has a 


(power) gain defined as 



g 




DirBction of wavG 
^^T^ropagation 


w. 


Wave motion in ail planes 
(a) Non -polarized wave 



y Wave motion in vertical plane 

only 

{b) Plane -polarized wave 


A 



(c) Cross -polarized waves 


AA and BB are orthogonal planes 
in which wave motion takes place 


FIGURE 1.2 Polarization of signal carrying waves 


Transmitter Repeater 


Repeater Receiver 



FIGURE 1.3 Arrangement of repeaters in a (Cable) transmission system 




Since the possible values of g may be very lai'ge, it is 
convenient to express gain in logarithmic units i.e. 


9 


dB 


10 log^Q g . 


Sometimes the signal strength <P> itself is expressed 
in logarithmic units, such as. 


dBW 


10 log,„( 


10 1 W 


) 


Any passive transmission system has a power loss rather 
than gain, since *^out ^ ^in' preferred to work with 
transmission loss or at tenv.ation defined as. 


L 



P. 

1 n 


out 


and 


‘"dB ■ " ^dB ” 

out 

Signal transmission by cables is very lossy : a lot of 
repeaters (mid-link power amplifiers where the signal power along 
with the noise content gets amplified) are needed. (See Figure 
1.5) On the other hand radio propagation (Microwave is at the 
upper extreme of the Radiowave portion of the spectrum) reduces 
the required number of repeaters and eliminates the cost of 
cables. Furthermore, large bandwidths are available for use in 
the microwave range which means more information can be carried 


by the communication channels. 



A lins-o/sight propa.ga.iton means the radio wave 
travels a direct path from the transmitting to the receiving 
antenna. The path, loss or the fr&&-spac& loss on 1 ine-of -s i ght 
path is due to spherical dispersion of the radio wave. This loss 
is given by 


L 


A. 




where 


- xpress 1 ng 


X 

/ 

I 

c 

in km and / in GHz ( range 


= signal wavelength; 

= signal frequency; 

= path length; and 

= spe ed of light, 

of Microwave ), we get 


L 


dB 




DirBctional ant&nncis have a focussing effect that acts 
like amplification , t.e. 


P 


out 


P , where 

L in 


and Qp are antenna gains at the transmitter and receiver 
respectively and L is the path ioss. The maximum transmitting or 
receiving gain of an antenna with effective aperture area is 

4nA 

Q = , where 

c = 5 X lo' Km/s. 

he value for A for a horn or dish antenna approximately equals 
t's physical area, and large parabolic dishes may provide gains 
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in excess of 60 dB. As described under Section 1, the focussing 
effect of the directional antenna also innplies attenuation in 
directions other than the antenna direction. 

3 Interference Mechanisms 

Interference refers to the contami nat i on of an 
information bearing signal by another similar signal, usually 
From a human source. This occurs in radio communication when the 
'eceiving antenna picks up two or more signals in the same 
'requency band. Interference may also result from multi -path 
iropagation. Regardless of the cause, severe interference 
prevents successful recovery of the message information. 

Propagation model of interference 

The interference model is built on a representation as 
hown in Figure 1.4. The RF-decoup 1 i ng is calculated by the 
ormula below 

= D + + ap + AG Aa + <5 . ( I . 1 > 


nere 


0 


ZO log 


D . 


0 


Xi>S 


difference between 


propagation loss in interfering and wanted signal 


paths ; 



Ill 



G = antenna gain 

a = attenuation between 
antenna and transmitter output 

D = path length 


D 

D 

D 


u , Daj : Wanted signal paths 


lA ’ 

2A ’ ^ 2.B 


Interfering signal paths 


A,B ; MW stations of interfering 

hop 

1,2 : MW stations of sample hop 


FIGURE L4 Propagation model of interference 
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- angle attenuation at interfering transmitter; 

= angle attenuation at interfered receiver; 

= difference between antenna gains at wanted and 
interfering signal transmitters; 

= difference between attenuations at wanted and 

interfering signal transmitters; 

= additional obstructions in interfering signal 

path . 

The above model of interference is based on the 
propagation formulae of signal waves. It does not characterize 
•.he type of interference taking place and also does not know 
whether a particular type of interference will actually take 
Place. Such a model can be referred to as the pTopa.ga.ttoT>. model 
>/ inter ference . 

There are three types of interference possible i.e'. 

(i) Co-channel interference 
(•11) Adjacent-channel interference 

( 1 i 1 )Rece i ver/Transmi tter Intermoduiat i on interference, 
nterference calculations, based on the propagation model, must 
herefore be made for each of the above possibilities. 

. 2 Co-channel model of" int-erference 

Both the desired and the undesired signals are carried 
n the same frequency channel and the receiver is unable to 


T 

AG 


Aa 



A = that atnount in dB by which the side-band noise is below 
the carrier power; 

B and C are constants; and 

AF = frequency offset from the carrier wave frequency at 
which the sideband noise strength is to be estimated. 
Such an equation has bo be constructed for each carrier wave 
frequency range. Obviously, the value for A depends upon the 
efficacy of filtering carried out at the transmitter. 
Adjacent-channel interference exists if the value of A is smaller 
than a specified threshold value. The extent of ad jacent -channe 1 
interference can also be calculated by calculating the amount of 
interference based on the the propagation model and then applying 
to it the Interference Reduction Factor <IRF) which depends upon 
the frequency offset <AF>. 

3. -4 Rocwivor /Transmit. tor Intormodulation intorf oronc© model 

This can occur at both the transmitter and the receiver 
and is a result of nonlinear behaviour of the circuitry at the 
transmitter output/rece i ver input stages. Interactions occur 

oetween the various signals present at a site because the 
circuitry is responsive to all of these signals, though not 
iqua 1 1 y so . 

Receiver i nt e rmodu 1 at i on occurs when a receiver 
■requency is expressible as an integer linear sum of the 
nterfering signal frequencies ii.e. coefficients of the linear 
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terms in the sum are all integers). The sum is also called the 
intermodulat ion product frequency and the sum of the absolute 
values of the integer coefficients is called the order of the 
intermodulation product. The amplitude of the product wave 
depends upon it’s order as well as it’s frequency offset from the 
victim receiver’s frequency. Increased linearity in the 
circuitry at the receiver input stage helps in reducing the 
intermodulation interference. 



Transmitter intermodulat ion 

occurs 

when 

a 

s ignal 

rece i ved 

from one 

transmi tter 

gets 

comb i ned 

wi th 

the 

output 

signal of 

a second 

transmi tter 

in it's 

output 

stage 

to 

produce 


1 ntermodu lat i on product frequencies which are radiated by the 
second transmitter. Selective f i 1 ters/ i solators at the output 
stage of the second transmitter is useful for reducing the 
amplitude of the signal received from the first transmitter. As 
the path loss of the involved signals is double that for a 
receiver i nt ermodulat i on , the transmitter intermodulation 
products are considerably weaker in amplitude. 

Intermodulation interference is said to occur if the 
amplitude of the i ntermodu lat i on product is greater than the 
specified threshold value at the victim receiver. 

4 Intornational>^National Roguiations for Frequency Usage 

International agencies like CCIR and ITU ( Internat lonal 
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Telecommunication Union) whose membership comprises of almost all 
the nations of the world, administer the usage of electromagnetic 
spectrum. The world has been divided into various regions for 
the purposes of regulating frequency usage. Each region has an 
international allocation of frequency bands/channe 1 s for use in a 
specified type of telecommunication service. The national 
allocations are also made similarly for the distribution of 
spectrum resources amongst the various users. The frequency plan 
prepared by each user (for *his' network) must be registered with 
the National Frequency Registration Board for proper 
documentation and administration of the national frequency 
a 1 1 ocat 1 ons . 

Keeping in view the above issues the frequency planner 
should be aware of all the relevant regulations governing 
frequency usage in 'his* geographical zone and also in the type 
of telecommunication service ’his’ networks intends to provide. 



APPENDIX II 


1 The Indian Railways Microwave Network 

The Indian Railways has got it's own microwave network, 
spanning the entire nation, to meet it's communication needs. 
There are about 100 base stations < transmi tter/rece i ver sites) 
and about 500 hops (communication links) spread almost uniformly 
ail over India. Presently, the network is being modernized from 
operating analog systems to using digital microwave systems. 

As per the plan prepared by the Indian Railways, there 
will be three distinct types of networks during three different 
time per 1 ods : - 

(i) Existing Analog Network - The communication links are 
designed to handle only analog signals. This network will be in 
use till the implementation of the digital microwave systems. 

( 11 ) Planned Digital Network - New digital systems will be 
deployed along with the retained analog systems for the microwave 
1 1 nks . 

<iii)Future Digital Network - This will be predominantly 
digital with some analog systems retained. 

From the above description it is clear that there is a 



transition from using analog systems to operating digital systems 
and the transitional phase will witness analog and digital 
systems co-existing in the planned digital network . This is 
because the analog systems cannot be shut down abruptly to deploy 
the digital systems - care has to be taken that the links ‘are not 
shutdown, as the telecommunication needs of the Indian Railways 
ajre of a continuous nature. 

2 Objectives of Frequency Planning 

The objectives pursued by Indian Railways for the 
frequency planning of it’s microwave network are long-term in 
nature and they are broadly three-fold :- 

<i) There should be a stable, dominantly digital microwave 
network to meet the Indian Railways telecommunication 
needs . 

(ii) The transition from the existing network operations to 
the future network operations must be smooth, 
especially during the interim planned digital network 
phase when digital and analog systems must co-exist. 
The frequency plan for the planned network should be 
such that little alteration is required in the future 
network and the shutdown cases for the analog systems 
in the planned digital network are kept to the minimum. 



< i i i ) Frequency planning for the planned digital network is 
most crucial and must be carried out af,ter taking into 
consideration the requirements of the future network. 

3 Process chart- for Manual Frequency Planning 

The process chart in Figure II.l illustrates the 
methodology followed for manual frequency planning. The Survey 
Map shows all the base stations with their coordinates <i.e. 
latitude and longitude) while the Sample Hop Data pertains to the 
technical details of the hop under study for interference from 
other hops within i t.* s coordination area. The Concept is based 
upon the formulae for compatibility calculations and it also 
gives guidelines for frequency allocations based on mandatory 
regulations and experience. It also stipulates the minimum 

requirement of RF-decoup 1 i ng for each link in order to avoid 
interference from other links. 

The data of the sample hop is marked on the survey map 
for ready reference. A Route Map is drawn, showing the various 
interfering hops (within the coordination area around the sample 
hop) for which compatibility must be established. 

The frequency arrangements are made for the sample hop 
with the polarization direction also being determined. The 

frequency channels are chosen such that compatibility 
calculations are most likely to succeed. 




FIGURE 11.1 Manual Process of frequency planning 
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An Engineering table is prepared for ready reference 
while making compatibility calculations. Besides the system 
operation data, angles of RF-decoupl i ng and lengths of the wanted 
and the interfering signal paths are also included in it. 


For ease of 

evaluating the RF-decoupl ing 

value 

<^RF> 

an 

aux i 1 iary tabl e 

i s 

prepared using which 

the 

^RF 

value 

i s 

calculated and 

the 

results are evaluated 

to 

check 

if 

the 


compatibility is satisfied as per the concept. 

If the compatibility is not established then the 
iteration is repeated afresh by changing the frequency 
arrangements. Once compatibility is secured, the frequency 
arrangement is documented and registered with the National 
Frequency Registration Board for the purposes of information and 
administration. 

4 Computerized Frequency Planning Process 

The computerization of the frequency planning process, 
as outlined in the Indian Railways report C*3 , is to be carried 
out using a mainframe computer. It is for most parts an 
automation of the existing manual process of frequency planning. 
It involves creation and maintenance of large databases for the 
system models, the base stations, and the hops. The various 
fields of data stored in each of these databases are as given 



below 


(i) Model - Microwave System used (FM 120/FM 500/FM 

960/ORS 54) 

“ Antenna Type (Standard Performance/High 
Performance) 

- Compatibility requirements 

<ii) Station - Altitude of site 

- Tower height 

- Coordinates of site 

( i i i )Hop - System 

- Antenna height 

- Antenna feeder attenuation 

- Power attenuation 

- Frequency arrangement 

- Polarization 

- Number of frequency channels required 

- Distinct network (Existing/Planned/Future) 

- Data of frequency channels 

Some of the field values must be filled in two records 
for each hop - one for each end station. 

System operation data is available for each microwave, 
system including information regarding each of the following 

- Co-channel decoupling required 

- Minimum angular separation required for each 
antenna type 

- Polarization permissible 
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- Frequency range and spacing 

- Attenuation at various stages in the link 

- IRF curves 

The compatibility calculations are made by a computer 
program called the HOP-ENGINEERING routine. It is called by 
specifying the Hop keys (two keys - one for' each station) and the 
System key. The calculations are made using a compatibility 
natrix which is already calculated and stored. The output is a 
listing of those critical cases where compatibility is not 
satisfied. The database updation is carried out using a separate 
:omputer program called the HOP-DIALOGUE routine. 



Discussion during the Thesis-Examination 


1. Since the intermodulation evaluation module consumes a lot 
of time in iterations and it’s effectiveness is also dependent 
upon the skills of the user, it was felt that subsets of 
1 ntermodulat i on-free frequencies should be generated prior to 
entering the coloring module by using the same nonlinear 
programming technique. Alternatively, a heuristic suggested was 
to group frequencies such that the ratios between any two 
frequencies are certain types of numbers (for example, 0.1?, 
0.8?, 0.7 etc.). This might result in the group of frequencies 
to be i nt ermodu 1 at i on- f re e . 

2. An interesting observation was that instead of partitioning 
the graph GC to be colored using the property of minimum cut, one 
could try to see if some "n i ce 1 y-co lorab 1 e" subgraphs can be 
obtained by deleting some vertices of GC . For instance, if 
planar subgraphs are obtained, they are easily ?-colorable (in 
principle, they are 4 -co 1 orabl e ) . This property would be 
desirable if the subsets of frequencies which are 
i n t e rmodu 1 at i on - f r e e are of relatively small cardinality (? to 
10, for example). Once the subgraphs are colored, merging occurs 
by re- 1 nt roduc i ng the vertices deleted earlier with a suitable 
coloring, one-by-one. 

The methodology of recognizing "n i ce 1 y -co lorab 1 e" 

subgraphs has to b© exa-rinned fui'^ther. 
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